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ABSTRACT 


The  results  of  theoretieal  research  on  Free  Electron  Lasers  (FELs)  are  presented. 
Basic  FEE  physics  is  reviewed,  using  a  previously  developed  classical  theory.  Numerical 
simulations  based  on  this  theory  are  described,  and  numerous  examples  show  how  they 
have  been  used  to  increase  understanding  of  existing  FELs  and  to  help  plan  new 
experiments.  Single-mode  simulations  that  follow  the  evolution  of  a  single-frequency 
plane  wave  provide  insight  into  important  physical  effects  in  FELs.  Results  show  how 
these  simulations  are  used  to  evaluate  new  FEE  designs  such  as  inverse-tapered  and  step- 
tapered  undulators.  Longitudinal  multimode  simulations  model  plane  waves  using  finite- 
length  electron  and  optical  pulses.  These  simulations  are  used  to  study  coherence 
evolution  in  various  FEE  designs,  and  to  explain  effects  such  as  limit-cycle  behavior. 
Transverse  multimode  simulations  that  allow  for  the  finite  transverse  dimensions  of  the 
optical  wavefronts  include  the  effects  of  optical  mode  distortion.  These  simulations  are 
currently  being  used  to  design  short  Rayleigh  length  optical  cavities  that  are  less  sensitive 
to  mirror  damage.  Four-dimensional  simulations  are  also  described,  which  follow  the 
optical  wavefront  in  x,  y,  z,  and  t,  including  the  effects  of  multiple  longitudinal  and 
transverse  modes.  These  simulations  are  computationally  intensive,  but  may  play  an 
important  role  in  the  design  of  future  high-power  FELs. 
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LIST  OF  SYMBOLS 


Throughout  this  dissertation,  cgs  units  are  used,  except  for  the  dimensionless 
variables  mentioned  below,  which  are  defined  in  the  text. 

A  magnetic  vector  potential 

a  dimensionless  optical  field 

B  undulator  magnetic  field 

optical  magnetic  field 

P  dimensionless  relativistic  electron  velocity 
c  speed  of  light  in  vacuum 
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e  electron  charge 
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electric  field  of  optical  wave 
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fiy)  electron  phase  velocity  distribution 

(j)  phase  of  optical  wave 


G  optical  field  single-pass  gain 

Y  Lorentz  faetor 

T]  single-pass  extraction  efficiency 

J  current  density 

j  dimensionless  eurrent  density 

k  optical  wavenumber 

undulator  wavenumber 
L  undulator  length 

4  optieal  pulse  length 

A,  optical  wavelength 

A-o  undulator  wavelength 

m  electron  mass 

N  number  of  undulator  periods 

n  number  of  passes  through  the  undulator 

V  dimensionless  eleetron  phase  veloeity 

dimensionless  synehrotron  frequency 

P  dimensionless  optieal  power 

Q  resonator  quality  factor 

q  electron  phase  veloeity  shift  from  resonance 

Te  eleetron  beam  radius 

rm  dimensionless  mirror  radius 

Tc  dimensionless  radius  of  eurvature 

p  eleetron  density 
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S  resonator  cavity  length 

G  g  electron  phase  velocity  spread  due  to  Gaussian  energy  spread 

Gg  electron  phase  velocity  spread  due  to  angular  spread 

G^  dimensionless  electron  pulse  length 

t  time 

T  dimensionless  time 
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EXECUTIVE  SUMMARY 


The  free  electron  laser  (PEL)  is  of  interest  to  the  Navy  as  a  scientific  tool  and  a 
shipboard  weapon  system  to  destroy  sea-skimming  missiles.  The  FEE  uses  a  relativistic 
electron  beam,  which  passes  through  a  magnetic  undulator  inside  an  optical  cavity.  The 
electron  beam  must  be  high-energy  to  produce  short  optical  wavelengths.  The  magnetic 
undulator,  which  consists  of  spatially  alternating  magnetic  fields,  causes  the  relativistic 
electrons  to  “wiggle”  and  thus  radiate  and  amplify  co-propagating  light.  The  optical 
cavity  stores  some  of  the  light,  allowing  it  to  develop  a  narrow  spectrum  through  mode 
competition  over  many  passes. 

The  key  advantages  of  FEEs  over  existing  light  sources  are  tunability,  high 
power,  and  high  efficiency.  While  the  wavelengths  of  conventional  lasers  are  constrained 
by  the  allowed  energy  transitions  of  electrons  bound  in  atoms,  the  wavelength  of  an  FEE 
can  be  tuned  over  a  broad  range,  from  microwaves  to  X-rays,  since  the  electrons  are 
“free”  to  make  energy  transitions  over  a  continuum  of  states.  An  FEE  can  amplify  large 
optical  power  without  damage,  since  there  is  no  heat-absorbing  medium  within  the 
optical  cavity.  FEEs  can  be  highly  efficient:  up  to  40%  extraction  efficiency  from  the 
electron  beam  has  been  demonstrated,  and  up  to  10%  wall-plug  efficiency  may  be 
possible. 

The  first  part  of  the  dissertation  explains  the  classical  theory  of  FEEs.  The 
resonance  condition  relates  the  optical  wavelength  to  the  electron  beam  energy  and 
properties  of  the  magnetic  undulator.  The  FEE  pendulum  equation  describes  the 
evolution  of  the  electrons  with  respect  to  the  optical  and  undulator  fields.  Phase-space 
plots  are  used  to  show  how  electrons  bunch  on  the  scale  of  an  optical  wavelength, 
enabling  power  exchange  between  the  electrons  and  the  optical  field.  The  FEE  wave 
equation  describes  the  evolution  of  the  optical  field  due  to  electron  bunching.  Examples 
show  how  the  electron  beam  quality  (energy  spread  and  emittance)  affect  bunching  and 
FEE  gain.  As  strong  optical  fields  develop,  the  electrons  begin  to  undergo  synchrotron 
oscillations  in  phase  space,  leading  to  saturation  of  the  optical  power.  Tapering  the 
magnetic  field  along  the  undulator  length  can  extend  the  saturation  limit. 
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The  following  chapters  describe  the  results  from  a  set  of  numerical  PEL 
simulations  based  on  this  classical  theory.  All  of  the  simulations  use  the  pendulum 
equation  to  follow  the  electron  evolution,  and  the  wave  equation  to  follow  the  optical 
field  evolution.  These  simulations  are  described  in  order  of  increasing  complexity.  The 
simplest  ones,  single-mode  simulations,  model  the  interaction  of  a  single-frequency  plane 
wave  with  the  relativistic  electron  beam  as  it  travels  through  the  undulator.  These 
simulations  provide  insight  into  important  physical  effects  in  FELs,  and  are  used  to 
evaluate  new  FEE  designs  such  as  inverse-tapered  and  step-tapered  undulators. 

Eongitudinal  multimode  simulations  use  finite-length  electron  and  optical  pulses 
to  model  plane  waves  of  multiple  frequencies  inside  the  resonator.  The  slightly  slower- 
moving  electron  pulse  slips  back  relative  to  the  optical  pulse,  reducing  the  overlap  and 
the  optical  power  gain.  These  simulations  are  used  to  study  coherence  evolution  in 
various  FEE  designs,  and  to  explain  effects  such  as  limit-cycle  behavior. 

Transverse  multimode  simulations  allow  for  the  finite  transverse  dimensions  of 
the  optical  wavefronts.  The  parabolic  wave  equation  is  used  to  include  the  effects  of 
optical  diffraction.  Diffraction  generally  reduces  the  overlap  between  the  optical 
wavefront  and  the  electron  beam,  but  under  special  conditions  the  FEE  mechanism  can 
continuously  distort  the  optical  wave,  improving  the  overlap  and  enhancing  the  gain. 
These  simulations  are  useful  in  studying  new  optical  cavity  designs. 

Four-dimensional  simulations  follow  the  evolution  of  the  electrons  and  the  optical 
wavefront  in  x,  y,  z,  and  t,  modeling  multiple  longitudinal  and  transverse  modes.  In 
general,  both  longitudinal  and  transverse  modes  are  present  and  interact  with  each  other. 
These  simulations  include  the  effects  of  longitudinal  pulse  slippage  and  transverse  mode 
distortion.  However,  they  are  computationally  intensive,  requiring  large  amounts  of 
memory  and  computer  time.  Four-dimensional  simulations  may  play  an  important  role  in 
the  design  of  future  high-power  FEEs,  which  require  short  Rayleigh-length  optical 
cavities  that  are  less  sensitive  to  mirror  damage. 
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I.  INTRODUCTION 


This  dissertation  presents  the  results  of  Free  Eleetron  Laser  (FEE)  research  over 
the  past  10  years  at  the  Naval  Postgraduate  School  (NPS).  The  author  has  been  deeply 
involved  in  all  of  this  work,  making  significant  contributions  as  the  co-author  of  10 
articles  published  in  refereed  journals,  and  as  the  lead  author  of  7  of  those  articles. 

An  FEE  is  a  new  type  of  laser,  which  uses  a  relativistic  electron  beam  to  create 
and  amplify  coherent  radiation.  In  a  conventional  laser,  the  electrons  in  the  lasing 
medium  are  bound  in  atoms,  and  a  particular  transition  between  fixed  atomic  energy 
levels  determines  the  radiation  wavelength.  In  an  FEE,  the  electrons  are  “free”  to  make 
an  energy  transition  over  a  continuum  of  states,  producing  radiation  over  a  broad  range  of 
wavelengths  from  microwaves  to  X-rays,  depending  on  the  design  parameters.  Since 
there  is  no  heat-absorbing  medium  inside  the  optical  cavity,  an  FEE  can  process  large 
amounts  of  optical  power  without  damage.  At  NPS,  we  do  theoretical  research  and 
numerical  modeling  of  FELs,  helping  to  design  new  FELs  at  many  universities  and 
laboratories  around  the  world,  and  to  improve  understanding  of  existing  experiments. 

Chapter  II  of  this  dissertation  reviews  the  basic  physics  of  FELs  using  a 
previously  developed  theory  [1],  explaining  the  FEE  interaction  in  terms  of  the  classical 
pendulum  equation  and  Maxwell’s  wave  equation.  The  equations  are  expressed  using 
dimensionless  parameters,  which  give  insight  into  the  FEE  mechanism,  and  facilitate 
numerical  modeling.  We  have  used  this  theory  to  develop  simulations  of  FELs  on 
desktop  computers.  These  simulations  have  been  benchmarked  against  theoretical 
formulas  and  by  comparison  to  numerous  experiments.  They  have  also  provided  the 
basis  for  FEE  simulations  used  by  other  researchers. 

The  chapters  that  follow  present  the  various  types  of  simulations  we  have 
developed;  single-mode,  longitudinal  multimode,  transverse  multimode,  and  four¬ 
dimensional  simulations.  The  basic  structure  of  the  simulations  is  described,  including 
their  assumptions  and  limitations,  the  input  parameters,  and  the  output  they  produce. 
Examples  are  shown  of  how  we  have  used  these  simulations  at  NPS  to  model  FEE 
experiments,  to  study  important  physical  effects,  and  to  extend  the  basic  theory  to  new 
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FEL  concepts.  Current  NFS  research  is  highlighted,  including  our  ongoing  collaboration 
with  the  Thomas  Jefferson  National  Aecelerator  Faeility  (TJNAF)  to  develop  a  high- 
power  FEE  for  naval  applieations. 

Chapter  III  deseribes  our  single-mode  simulations,  whieh  model  a  single¬ 
frequency  plane  wave  interaeting  with  many  eleetrons  in  the  relativistic  electron  beam. 
We  have  used  these  simulations  to  study  the  effects  of  electron  beam  quality  on  FEE 
gain,  and  to  understand  basie  physieal  effects  in  new  FEE  designs.  This  research  helped 
to  design  the  SEEENE  FEE,  a  ground-based  laser  proposed  in  the  early  1990’s  for 
beaming  power  to  satellites.  In  recent  years  we  have  used  single-mode  simulations  to 
study  new  designs  for  the  current  and  proposed  FEEs  at  TJNAF,  resulting  in  several 
conference  presentations  and  a  published  artiele. 

Chapter  IV  describes  our  longitudinal  multimode  simulations,  whieh  model  plane 
waves  using  finite-length  eleetron  and  optical  pulses.  These  simulations  have  helped  to 
understand  new  effeets  such  as  limit-cycle  behavior  in  short  pulse  FEEs.  We  used  these 
simulations  in  the  early  1990’s  to  model  eleetron  beam  energy  modulation  in  the  SCA 
FEE  at  Stanford.  Our  longitudinal  multimode  simulations  contributed  to  the  design  of  a 
high-power  FEE  proposed  at  Boeing  in  the  mid- 1990’s,  predicting  the  development  of  a 
ehaotic  optical  spectrum.  More  reeently  we  have  used  these  simulations  to  study  mode 
competition  and  coherenee  evolution  in  the  proposed  X-ray  FEE  at  SEAC.  They 
eontinue  to  play  an  important  role  in  our  FEE  design  studies  for  TJNAF.  The  results  of 
this  work  have  been  presented  at  several  conferences  and  in  three  published  articles. 

Chapter  V  deseribes  our  transverse  multimode  simulations,  whieh  allow  for 
optieal  wavefronts  of  finite  transverse  dimensions.  We  have  used  these  simulations  to 
study  optical  mode  distortion  in  the  SEEENE  and  Boeing  FEEs,  and  they  have 
eontributed  to  the  design  of  new  optical  cavities  for  the  TJNAF  FEE.  These  researeh 
results  have  been  presented  at  several  conferences  and  in  two  published  artieles. 

Chapter  VI  deseribes  our  four-dimensional  simulations,  which  model  the  electron 
beam  and  the  optical  wavefront  in  x,  y,  z,  and  t,  ineluding  the  effeets  of  multiple 
longitudinal  and  transverse  modes.  These  simulations  are  computationally  intensive, 
requiring  large  amounts  of  memory  and  computer  time.  In  the  early  and  mid- 1990’s  we 
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developed  versions  of  these  simulations  that  eould  run  on  desktop  workstations,  but 
applied  only  to  FELs  with  narrow  electron  beams.  We  used  these  simulations  to  study 
coherence  evolution  and  optical  mode  distortion  effects  in  the  Boeing  FEE,  and  to  help 
design  an  infrared  FEE  at  TJNAF.  We  presented  results  of  this  research  at  several 
conferences,  and  in  three  published  articles.  We  are  currently  studying  new  numerical 
techniques  to  reduce  the  computational  requirements  for  our  four-dimensional 
simulations,  which  may  be  useful  in  the  design  of  high-power  FEEs  for  ship  defense. 
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II.  BASIC  FEL  THEORY 


A  Free  Electron  Laser  (FEL)  is  a  device  used  to  create  powerful,  tunable, 
coherent  radiation  for  scientific,  medical,  industrial  and  military  applications.  Figure  1 
shows  a  schematic  of  an  FEL  oscillator,  which  was  first  described  by  John  Madey  in  the 
early  1970’s  [2]. 


electron  beam 


resonator  mirror 

Figure  1.  FEL  Oscillator 


The  major  components  are  a  relativistic  electron  beam,  a  magnetic  undulator,  and  a 
resonator  cavity.  A  linear  accelerator  or  a  storage  ring  provides  the  relativistic  electron 
beam,  with  typical  energy  of  100  MeV  and  peak  current  of  100  A.  The  undulator  is 
either  a  series  of  alternating  permanent  magnets  or  a  helically  wound  electromagnet,  with 
typically  50  periods,  each  a  few  centimeters  long.  The  undulator  field  strength  is  on  the 
order  of  a  few  kilogauss.  The  resonator  consists  of  an  evacuated  optical  cavity  typically 
about  10  meters  long,  terminated  by  mirrors  on  each  end.  One  of  the  mirrors  is  partially 
transmitting  to  outcouple  the  light  for  applications.  The  electron  beam,  after  passing 
once  through  the  undulator,  can  either  be  directed  to  a  beam  dump,  or  recirculated  back 
to  the  accelerator  to  enhance  overall  efficiency  and  reduce  the  energy  of  the  dumped 
beam. 

In  an  FEL  amplifier,  there  are  no  mirrors.  Radiation  from  an  external  source  is 
amplified  over  a  single  pass.  It  is  also  possible  to  create  and  amplify  light  over  a  single 
pass,  with  no  external  source,  in  a  process  known  as  self-amplification  of  spontaneous 
emission  (SASE). 
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The  magnetic  undulator  gives  the  relativistic  electrons  a  slight  “wiggling”  motion 
from  side  to  side,  yielding  spontaneous  radiation  in  a  narrow  cone  in  the  forward 
direction.  The  radiation  wavelength  is  determined  by  the  electron  beam  energy,  the 
undulator  wavelength,  and  the  magnetic  field  strength.  The  transverse  motion  of  the 
electrons,  when  properly  coupled  to  the  transverse  electric  field  of  the  radiation,  produces 
stimulated  emission,  amplifying  the  light  as  it  passes  through  the  undulator. 

FELs  offer  significant  advantages  over  existing  radiation  sources.  Unlike  a 
conventional  laser,  the  available  wavelengths  are  not  restricted  by  allowed  transitions 
between  energy  levels  of  bound  electrons.  An  FEE  can  be  designed  to  operate  at  any 
wavelength  from  microwaves  to  X-rays,  and  a  particular  FEE  can  also  be  tunable  over  a 
broad  band.  A  wavelength  range  of  xlO  for  a  single  FEE  has  been  demonstrated  [3],  and 
a  range  of  xlOO  may  be  possible.  Unlike  a  synchrotron  light  source,  FEE  radiation  is 
highly  coherent.  Since  there  are  no  atoms  or  molecules  inside  the  optical  cavity,  only 
relativistic  free  electrons,  exhaust  “heaf  ’  is  removed  at  the  speed  of  light,  so  it  is  possible 
to  process  megawatts  of  optical  power  inside  the  cavity.  An  FEE  can  be  highly  efficient; 
up  to  40%  energy  extraction  from  the  electron  beam  has  been  demonstrated  [4].  FEEs 
can  also  produce  very  short  (picosecond)  optical  pulses,  which  are  useful  for  many 
applications. 

The  first  FEE  was  built  at  Stanford  in  the  mid  1970’s  [5].  Since  then,  hundreds 
have  been  constructed  around  the  world,  with  wavelengths  ranging  from  0.1  pm  to 
several  hundred  pm  [6].  The  most  powerful  FEE  is  at  Jefferson  Eabs  in  Virginia  [7];  it 
has  produced  2  kW  of  average  power  at  3  pm,  and  will  soon  be  upgraded  to  10  kW  [8]. 
There  are  proposals  to  build  X-ray  FEEs,  with  wavelengths  from  1  to  10  A,  at  DESY 
(Germany)  and  SLAC  (Stanford);  these  would  be  the  world’s  first  X-ray  lasers,  with 
many  novel  applications  [9]. 

A,  RESONANCE  CONDITION 

The  primary  factor  that  determines  the  operating  wavelength  of  an  FEE  is  the 
electron  beam  energy  y  me  ,  where  y  =  ( 1  -  P  1  is  the  Eorentz  factor,  v  =  P  c  is  the 
electron  velocity,  m  is  the  electron  rest  mass,  and  c  is  the  speed  of  light.  In  the  reference 
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frame  of  the  eleetron,  the  undulator  wavelength  is  Lorentz  contraeted  to  =Xq/j  , 
and  the  eleetron  will  oscillate  and  emit  light  at  that  same  wavelength,  X'  =  Xq.  In  the 
laboratory  frame,  the  light  is  Doppler  shifted  to  the  wavelength  X=X' jly  ,  which  gives 


(1) 


This  equation  is  approximate  because  it  neglects  the  transverse  motion  of  the  electron  (it 
assumes  the  undulator  field  is  very  small),  but  it  shows  the  relationship  between  electron 
energy  and  optical  wavelength.  Thus,  to  achieve  very  short  wavelengths,  such  as  the 
proposed  X-ray  PEL  at  SLAC,  requires  a  high-energy  electron  beam  (>10  GeV,  so  that 
Y  >  20000),  whereas  infrared  radiation  can  be  produced  with  a  more  modest  accelerator 
( ~  50  MeV,  Y  ~  100  ).  Visible  and  ultraviolet  FELs  require  beam  energies  on  the  order  of 
500  MeV(Y  -1000). 


A  more  precise  relationship  can  be  derived  by  considering  a  single  electron 
traveling  at  speed  P^c  along  the  undulator  (z)  axis,  with  P^  « 1 ,  and  a  co-propagating 
photon  traveling  at  speed  c.  The  undulator  gives  transverse  motion  to  the  electron, 
enabling  energy  exchange  through  the  relativistic  Eorentz  equation, 

^  =  (2) 

at  me 


where  e  =  e  is  the  electron  charge  magnitude  and  E  is  the  electric  field  of  the  light 


wave.  Since  the  light  wave  is  transverse, E^=0  and  only  the  transverse  components  of 
the  velocity,  Pj^  =  P^^,  contribute  to  energy  exchange  in  (2).  For  optimum  energy 


exchange,  these  components  should  oscillate  at  the  same  frequency  as  E .  This  leads  to 
the  following  resonance  condition;  one  wavelength  of  light  passes  over  an  electron  as  the 
electron  travels  through  one  undulator  period.  In  the  corresponding  time  interval  At,  the 
electron  travels  a  distance  X^  =cP^At  and  the  photon  travels  a  distance  X^+X  =cAt. 


Thus, 


^1-p.  ^ 


V  Hz  y 


A-o . 


(3) 
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If  we  assume  a  highly  relativistie  eleetron  beam  (y  »  l)  and  ignore  the  transverse 
motion  of  the  eleetron  (P^  «cl),  then  P^  «1-1/2y^,  and  the  resonanee  eondition  (3) 
reduees  to  the  previous  result  (1). 


B,  ELECTRON  DYNAMICS 
1.  Betatron  Motion 

To  understand  the  proeesses  of  spontaneous  and  stimulated  emission  in  an  PEL,  it 
is  necessary  to  study  the  motion  of  the  electrons  in  the  undulator.  First,  we  will  derive 
the  electron  trajectories  in  a  linearly  polarized  undulator  with  no  light  present.  The 
magnetic  field  for  an  ideal  linear  undulator  can  be  written  in  rectangular  coordinates  [10] 

=.S(0,  sin (A:oz) cosh (A:oy),  cos(A:oz)sinh(A:oy))  (4) 


where  =27t:/A,q  is  the  undulator  wavenumber,  and  is  the  undulator  period.  Along 

the  undulator  (z)  axis,  the  field  is  given  by  Bsm{k(^z)y  with  an  rms  value  =  b/42  . 

Typically  kQy  «  1 ,  since  ~  5  cm  and  the  beam  radius  y  ~  1  mm.  Off  axis,  the  average 

field  B  =  B^{l  +  kQ^y^  -\ — j  increases  withy,  focusing  the  electrons  back  towards  the 

axis.  Since  the  field  (4)  is  independent  of  x,  there  is  no  natural  focusing  in  the  x 
direction.  In  practice,  the  beam  may  be  focused  using  parabolic  pole  faces  or  a 
quadrupole  lens. 


The  electron  motion  is  determined  by  the  relativistic  Lorentz  force  equations. 


4y) 

dt 


e 

me 


(E+px5) 


(5) 


plus  the  energy  equation  (2).  Using  (4)  for  the  undulator  field  and  assuming  no  light 
present  =  0  ,  y  =  oj ,  we  obtain 
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cB 

(3 [  cos  ( V )  j  sinh  ( ^0  J )  -  z  sin  (  )  cosh  ( ^0  j)  ] 

p  =  +  -^xcos(V)sinh(^oj) 
y  me 

= - i  sin  ( A:oz)  cosh  (A:oj) 

jmc 


(6) 


The  first  equation  is  easily  integrated: 


P^  =-V2 — eos(A:oz)cosh(A:oy) 


(V) 


where  K  =  eB^Xl^^  mc^  is  the  “undulator  parameter”,  a  dimensionless  quantity 

approximately  equal  to  1  for  most  FELs.  We  have  assumed  perfect  injection  in  the  x 
direction,  so  the  integration  constant  in  (7)  is  zero.  Now  substitute  (7)  into  the  last  two 
equations  in  (6)  to  get 

cK^ka 


Pv=- 


''0  2 

COS 


y 


(iCoz)sinh(2iCoy) 


cK^kn 


(8) 


p,  =  — ^sin(2V)cosh"(^oT)  • 


Averaging  over  many  undulator  periods  and  assuming  the  electron  is  near  the  undulator 
axis  [k^y^-l)  gives  ~ -{cK^kljy^^y  and  P^  «  0  .  Thus  the  averaged  equation  of 
motion  fory  is 

d^y 


^  Kk,L'^^ 


d\‘ 


\  Y  ) 


y 


(9) 


where  we  have  used  the  dimensionless  time  t  =  ct/T .  This  corresponds  to  simple 
harmonic  motion  at  the  dimensionless  betatron  frequency,  cOp  =  Kk^Ljy  =  InNKly  .  For 
typical  FEE  parameters,  cOp  «  27t:  ,  so  there  is  one  betatron  oscillation  along  the  undulator. 


The  general  solution  to  the  harmonic  oscillator  equation  (9)  ean  be  written 
y{x )  =  Jo  eos  (co  pT  j  + - sm  (co  pX  j 

CO  D 


(10) 
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where  and  0  are  the  initial  position  and  angle  of  an  eleetron  at  the  undulator  entranee 
(t  =  0)  ,  with  respect  to  the  undulator  axis.  The  initial  angle  is  related  to  the  initial  y- 

velocity  by  0  ,,  =>’(0)/l  where  (  )  indicates  a  derivative  with  respect  to  dimensionless 

time  T.  Equation  (10)  describes  the  electron  betatron  motion  in  a  linearly  polarized 
undulator  with  no  light  present. 

For  a  helical  undulator,  the  magnetic  field  can  be  written  in  cylindrical 
coordinates  using  Bessel  functions  [1].  Near  the  undulator  axis,  the  average  transverse 
field  can  be  expanded  for  to  get  =  b{\  + /4-{ — j,  which  increases  off 

axis  at  half  the  rate  of  the  linear  field,  so  the  betatron  frequency  is  reduced  by  a  factor  of 

V2  .  Since  the  field  is  axially  symmetric,  there  is  equal  betatron  focusing  in  both  the  x 
andy  directions. 

2,  FEL  Pendulum  Equation 

Now  we  consider  how  a  co-propagating  optical  wave  affects  the  electron 
dynamics.  For  a  helical  undulator,  assuming  that  betatron  focusing  keeps  the  electrons 
near  the  undulator  axis,  the  field  can  be  simplified  in  rectangular  coordinates  to 

Bjj  =  B [cos k^z, sin kQZ,0^ .  (11) 

Note  that  unlike  the  linear  case,  the  field  amplitude  is  constant  along  the  z-axis,  so 
B^^  =  B  in  the  definition  of  the  undulator  parameter  K. 

Now  introduce  a  helically  polarized,  monochromatic  plane  wave  with 
wavenumber  k  and  frequency  co,  traveling  in  the  +z  direction,  with  optical  fields  given  by 
Fj  =F'(cos\|/, -sin\|/,  0)  ,  =  F(sin\|/,  cos\|/,  O)  (12) 

where  \\f  =kz-cot  +  (j) ,  E  is  the  optical  field  amplitude  in  cgs  unit,  and  (j)  is  the  optical 
phase.  Integrating  (5)  using  the  undulator  field  (11)  and  the  optical  field  (12),  with  the 
assumption  y  » 1 ,  gives  the  electron  trajectories. 

Pi  =-— (coskoz,  sinkoz,  0)  ,  p,=l-^^^  .  (13) 

Y  2y 
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We  again  assume  perfect  injection  of  the  electrons  into  the  undulator,  so  the  constants  of 
integration  are  zero.  Imperfections  in  the  electron  beam  will  be  discussed  later.  Note 
that  the  amplitude  of  the  transverse  oscillations  is  small  since  Kjy  «cl,  but  as  we  will 
see  in  the  next  section,  it  is  these  oscillations  that  cause  spontaneous  emission,  and  couple 
the  electron  motion  to  the  optical  wave  to  enable  stimulated  emission. 


Using  the  expression  (13)  for  ,  the  resonance  condition  (3)  can  now  be  written 


X  = 


(14) 


Substituting  Pj^  from  (13)  into  the  energy  equation  (2)  gives 

dy  eKE 

Y  =  —  = - cos 

dt  y  me 


(15) 


where  ^  =(A:  +  A:o)z-cot  is  the  electron  phase  with  respect  to  the  combined  optical  and 
undulator  fields.  Equation  (15)  contains  the  physics  of  energy  exchange:  if  the  electron 
phases  C,  are  uniformly  distributed,  then  half  the  electrons  will  gain  energy  (y  >  0)  and 
half  the  electrons  will  lose  energy  (y  <  0) .  In  order  to  have  a  net  energy  transfer,  the 
electrons  must  be  bunched  on  the  scale  of  the  optical  wavelength. 


The  electron  phase  velocity  is  given  by  /Jt  =  (A:  +  )  P^c  -co  .  Taking  another 

time  derivative,  assuming  y  »  1  and  using  (13)  for  P^  gives 

dX 


dt 


-  kc 


ry] 

1  y'  J 

Iy  J 

(16) 


Substituting  y  from  (15)  and  using  the  resonance  condition  (14), 

dX  leKEX 


dt 


2 

y  m 


-cos 


ft-Kl-). 


(17) 


Assuming  the  electron  energy  y  mc^  doesn’t  change  much  over  an  undulator  period,  we 
can  simplify  this  equation  to  obtain 


»  dX  II  r;' 

V  =  |a|cos((,  +(})) 


(18) 
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where  v  =  Z[(A:  + A:o)  -  A:]  is  the  dimensionless  electron  phase  velocity,  and 
\a\  =  ATiNeKLEIy^mc^  is  the  dimensionless  optical  field  amplitude.  The  variable  v  is 

also  known  as  the  “resonance  parameter”,  since  v  =  0  for  an  electron  that  satisfies  the 
resonance  condition. 

The  result  (18)  also  holds  for  a  linear  undulator  if  we  replace  K  in  the  definition  of 
|a|  by  if  Jj  (^)) ,  where  andJj  are  Bessel  functions  of  the  first  kind,  and 

^  =K^  I2[\  +  K^).[\] 

Equation  (18)  is  known  as  the  FEE  pendulum  equation.  It  describes  the 
microscopic  motion  of  each  electron  along  the  z-axis  in  a  reference  frame  one  optical 
wavelength  long,  traveling  at  the  speed  of  a  resonant  electron.  In  that  reference  frame, 
the  electron  motion  is  analogous  to  a  classical  pendulum.  [11] 

To  visualize  this  motion,  it  is  useful  to  look  at  a  phase  space  plot  such  as  Figure  2 

below,  which  shows  the  electron  evolution  in  ^  and  v  for  a  typical  low-gain  FEE. 
This  figure  was  produced  by  numerically  solving  the  pendulum  equation  (18)  in  a  single¬ 
mode  simulation,  as  described  in  a  later  chapter.  Twenty  sample  electrons  are  injected  at 
resonance,  each  with  initial  phase  velocity  Vo=v(0)  =  0,  and  their  initial  phases  C,. 

evenly  distributed.  A  color  scale  is  used  to  show  their  evolution  from  the  beginning  of 
the  undulator  (r  =  0 ,  yellow),  to  the  end  of  the  undulator  (t  =  1 ,  red).  The  electrons  are 
seen  to  follow  pendulum  phase  space  orbits.  The  separatrix,  drawn  in  red  for  reference, 
is  defined  by  =  2|a|[l  +  sin(^^  +^)]  •  The  separatrix  passes  through  the  unstable  fixed 

points  at  ,v)  =  (-7t:/2,0)  and  (37t:/2,0) ,  with  peak-to-peak  height  4|a|’^^.  Inside  the 
separatrix,  the  electrons  follow  closed  paths  around  the  stable  fixed  point  at 
,v )  =  (7t:/2,0)  ,  corresponding  to  a  pendulum  swinging  back  and  forth;  outside  the 

separatrix  are  open  paths,  corresponding  to  a  pendulum  swinging  “over  the  top”.  Notice 
at  the  end  of  the  undulator  the  electrons  are  bunched  near  ^  =  7t:/2  . 
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Figure  2.  Electron  evolution  in  (C,v)  phase  space,  optical  field  gain  G(t) 
and  phase  (j)  (r )  for  a  single  pass  through  a  typical  low-gain  EEL, 
with  electrons  injected  at  resonance  (vg  =  O) .  The  electrons  follow 
pendulum  phase  space  orbits,  and  there  is  no  net  gain  (G  w  O) . 


A  change  in  phase  velocity  Av  can  be  related  to  an  energy  change  Ay  using  the 
definition  of  v  after  (18),  the  definition  of  (13),  and  the  resonance  condition  (14),  to 
obtain  Av  =  An  NAy  /  y  .  Thus  the  electrons  gain  energy  when  they  move  up  on  the 
phase  space  plot,  and  lose  energy  when  they  move  down.  In  figure  2,  half  of  the 
electrons  gain  energy  and  half  lose  energy.  The  net  result  is  very  little  energy  transfer,  as 
shown  by  the  gain  evolution  plot  G(t)  in  the  upper-right.  The  gain  is  defined  by 

G(t)  =  [p(t)-Po]/Eo,  where  P  (t)  is  the  optical  power  at  time  t  ,  and  Pq  =P{0)  is 
the  initial  optical  power.  The  lower-right  plot  shows  the  optical  phase  evolution  (j)  (t  ) , 
which  causes  the  separatrix  to  shift  slightly. 

figure  3  shows  another  simulation  with  the  electrons  injected  slightly  off 
resonance  at  Vg  ~2.6.  Some  electrons  now  follow  open  orbits  above  the  separatrix; 
those  on  the  left  side  of  the  plot  move  ahead  faster,  while  those  on  the  right  move  slower. 
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The  net  effect  is  bunching  near  C,  =tz  ,  where  the  electrons  “fall”  towards  smaller  values 

of  V,  losing  energy  to  the  optical  wave.  Thus  the  final  gain  increases  to  about  13%, 
whereas  in  the  previous  case  it  was  only  about  0.5%.  Also  note  that  the  optical  phase 
shift  has  decreased. 
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Figure  3. 
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Electron  evolution  in  ,v )  phase  space,  optical  field  gain  G  (t  ) 

and  phase  (j)  (r )  for  a  single  pass  through  a  typical  low-gain  EEL, 

with  electrons  injected  slightly  off  resonance  (vo=2.6).  The 

electrons  still  follow  pendulum  phase  space  paths,  but  now  there  is 
significant  energy  extraction  and  gain  of  about  13%. 


If  the  electrons  are  started  farther  off  resonance  (larger  Vq),  the  phase  space  orbits 
flatten  out  and  there  is  less  gain.  If  the  electrons  start  with  a  negative  phase  velocity 
(Vo<0),  they  bunch  near  C,  =0  and  absorb  energy  from  the  optical  field.  Eigure  4 
shows  the  gain  in  a  single  pass  through  the  undulator,  plotted  versus  initial  electron  phase 
velocity  Vg ,  generated  by  repeating  a  single -mode  simulation  while  varying  Vg .  The  gain 

curve  in  Eigure  4  is  antisymmetric,  and  the  peak  gain  is  about  13%  at  Vg  »  2.6 .  The  gain 
bandwidth  is  Av  «  27t:  .  Extracting  energy  in  excess  of  the  gain  bandwidth. 
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Av  =AtiN  Ay  ^  In  ,  reduces  the  gain,  so  the  maximum  theoretical  efficiency  of  an 

PEL  (single  pass  energy  extraction)  is  Ay  /y  ^IjlN  . 


-13% 


13% 


0.0 


Figure  4.  Single -pass  gain  versus  initial  electron  phase  velocity  for  a  low- 


gain  FEE  with  weak  optical  fields.  The  gain  curve  is  anti¬ 
symmetric,  with  peak  gain  of  about  13%  at  Vq  «  2.6 . 


Our  derivation  of  the  pendulum  equation  (18)  assumed  perfect  injection  of  the 
electron  beam  into  the  undulator.  If  an  electron  is  injected  off-axis,  at  a  displacement 

or  an  angle  0^, ,  it  will  undergo  betatron  motion  described  in  (10).  With  the  assumptions 
0^,  <K7t:  and  « 1 ,  it  can  be  shown  [12]  that  the  corresponding  phase  velocity  is 
reduced  by 


(19) 


Note  that  this  shift  is  constant  in  time,  so  with  an  adjustment  to  the  initial  phase  velocity, 
we  can  still  use  the  pendulum  equation  to  determine  the  electron’s  phase  space  evolution. 

In  a  realistic  beam,  the  electrons  enter  the  undulator  with  a  random  distribution  of 
positions  and  angles  characterized  by  rms  spreads  and  0^,  .  Each  of  those  can  be 

changed  using  focusing  fields,  but  their  product,  known  as  the  emittance  =yoQy  ,  is 

fixed  [1].  Experimentalists  try  to  “match”  these  quantities  so  that  the  beam  is  not  too 
large  and  does  not  expand  too  much  along  the  undulator.  When  the  matching  condition 

Kk^yg  xjQy  is  satisfied,  the  phase  velocity  spread  Avp  »  27i  AWA:8,,/y  is  minimized. 
Figure  4  shows  that  the  natural  gain  bandwidth  of  an  FEE  is  Av  ~  271  ,  so  for  optimum 
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gain  the  emittance  should  be  limited  to  8™“  ^jXjlnNK ,  which  corresponds  to  about  1 

mm-mrad  for  typical  PEL  parameters.  For  a  short  wavelength  FEE,  or  one  with  many 
undulator  periods  N,  the  required  emittance  can  be  significantly  less. 

C.  OPTICAL  FIELD  EVOLUTION 

1.  Spontaneous  Emission 

The  small  transverse  oscillations  (13)  of  the  relativistic  electrons  cause  them  to 
radiate  light  into  a  forward  cone  of  angular  width  on  the  order  of  y  '  [13].  If  A  <  1 ,  the 
radiation  cone  stays  in  a  detector  on  axis  at  infinity,  but  if  A  >  1  the  cone  periodically 
deflects  out  of  the  detector,  producing  higher  frequency  harmonics  [14].  This 
spontaneous  emission  is  required  for  start-up  of  an  FEE  oscillator  or  SASE  FEE. 
According  to  the  FEE  resonance  condition  (14),  the  fundamental  frequency  is 

CO  =  2y  j .  The  total  radiated  power  at  each  frequency  can  be  found  from  the 

relativistic  Earmor  formula,  which  reduces  to  a  simple  result  that  for  each  pass  through  a 
typical  undulator  (A  « l) ,  only  a  small  fraction  of  the  electrons  in  the  beam,  aK^  «  10^^ , 
contributes  a  photon  to  the  fundamental  mode,  where  a  is  the  fine  structure  constant  [1]. 


2.  Stimulated  Emission 

Now  we  consider  how  the  optical  wave  and  electrons  evolve  as  they  interact  in 
the  undulator.  First,  we  derive  a  simple  formula  for  the  gain  in  weak  optical  fields, 
|a|<sC7T:  .  We  assume  a  perfect  electron  beam,  with  all  the  electrons  at  the  same  initial 

phase  velocity  Vq  and  their  initial  phases  evenly  distributed.  If  we  assume  low  gain, 
a~aQ,  then  we  can  expand  ^  and  v  in  powers  of  the  initial  field  amplitude  : 


+v® 


(20) 


To  zeroth  order  in  (no  optical  field),  ^  =^q+VoT  and  =Vq.  To  first 


order  in  , 
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(21) 


C  {')  =  - ^ [cos(C 0  +V oT )  -  cos C 0  +V oT  Sin^  0  ] 

Vo 

=  — [sin(Co+VoT)-sinCo]. 

Vn 


Since  an  electron’s  energy  is  proportional  to  its  phase  velocity  v,  we  must  have 
(v)-Vo  0  in  order  to  have  a  net  energy  transfer,  where  (•••)  indieates  an  average  over 
all  the  electrons.  But  the  eleetrons  are  evenly  distributed  in  so  (21)  implies 

=  0  ,  and  thus  to  determine  the  gain  the  seeond-order  term  is  required  [1] 

2 

V  (")  =  -^  [- i ( cos(2C  0  +  2v  oT )  -  eos  2C  0 )  +  cos(v  qT  )  - 1  - V  qT  sin C  0  cos(C  0  +  V  qT  )]  .  (22) 
Vq 

The  first  term  inside  the  braces  averages  to  zero  over  all  the  electrons,  so  we  are  left  with 

(2)  ^  [cos(v oT  )  - 1  +  oT  sin(v  qT  )] .  (23) 

Vq 


Recalling  Av  =  47t:  Ay  /y  ,  we  can  write  the  average  energy  change  of  an  eleetron  as 


Ay  me 


2 


yme^  ((Av  )-Vq) 
AtiN 


Y^/  (2)\ 


(24) 


To  obtain  the  total  energy  change  of  all  the  electrons  in  a  volume  element  dV 
inside  the  laser  beam,  multiply  (24)  by  the  number  of  eleetrons  in  that  volume,  pFdV , 
where  p  is  the  electron  beam  density  and  F  is  the  “filling  factor”,  defined  as  the  eross- 
section  area  of  the  electron  beam  divided  by  the  cross-section  area  of  the  optieal  mode. 
From  energy  conservation,  all  of  the  energy  lost  by  the  eleetrons  goes  into  the  optieal 
field.  The  radiation  energy  in  dV  is  2E^  dVjSn  .  We  define  the  gain  of  the  optical  field 
as  the  fractional  increase  in  the  optieal  field  power  or  energy,  giving 


G  = 


(pFdV)(y  me  ^  /AnN'j(v^^^'j 
2E^  dViSn 


(25) 


Simplifying  and  using  (23),  we  obtain 

G  =  ^[2-2cos(VqT)-VqT  sin(vQT)l 
Vo 


(26) 
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where  j  =  %N{^eKKL)^  pF mc^  is  the  dimensionless  eurrent  density.  This  result  is 
valid  for  weak  fields  (a^  <K  tt  )  and  low  gain  (7  <  tt  ) .  Equation  (26)  deseribes  the  gain 
curve  shown  in  Figure  4,  which  was  obtained  by  numerical  simulation. 

3.  Coherence  Development 

The  electron  phase  velocity  v  is  a  function  of  the  optical  wavenumber  k,  so  that  a 
change  in  v  implies  a  change  in  the  frequency  of  the  optical  mode  inside  the  resonator. 
Thus  by  varying  v  and  calculating  the  EEL  response,  we  can  obtain  the  gain  or  power 
spectrum. 

The  spontaneous  emission  spectrum  ^(v)  is  symmetric  in  shape,  centered  at 
V  =  0  with  a  width  Av  «  Itt  .  [1]  In  weak  fields  and  low  gain,  each  mode  will  have  gain 
G(v)  given  by  (26),  antisymmetric  about  v  =  0  with  a  peak  at  v  =2.6.  Each  pass 
through  the  undulator,  the  optical  power  in  each  mode  Piy )  will  evolve  according  to 

AP(v  )  =  5(v  )  +  E(v  )  [G(v  )  - 1/0]  (27) 

where  Q  is  the  resonator  quality  factor;  for  0  »  1 ,  the  loss  on  each  pass  is  given  by  1/0  . 
Over  many  passes,  mode  competition  will  cause  the  linewidth  to  narrow,  as  shown  in 
Figure  5.  The  bottom  graph  shows  the  spontaneous  emission  spectrum  ^(v ) .  Above  that 

is  the  theoretical  gain  spectrum  G(v)  for  a  low-gain  FEE  (identical  to  Figure  4).  Above 
that  is  the  optical  power  spectrum  evolution  P{y,n)  n  =  100  passes,  showing  that  the 
initially  broad  spectrum  narrows  and  shifts  over  to  the  peak  of  the  gain  curve,  due  to 
mode  competition.  At  the  top  is  the  final  power  spectrum  (v ) .  Eventually,  the  optical 

linewidth  is  transform  limited  by  the  optical  pulse  length,  which  is  determined  by  the 
electron  pulse  length. 
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Figures.  Spontaneous  emission  speetrum  ^(v),  gain  speetrum  G(v), 

optical  power  spectrum  evolution  P{y,n),  and  final  power 

spectrum  (v )  for  a  typical  low-gain  FEL,  showing  coherence 

evolving  through  mode  competition  as  the  laser  line  narrows  over 
n  =  100  passes. 


4,  FEL  Wave  Equation 

A  more  general  formula  for  the  optical  wave  evolution  can  be  derived  by  starting 


with  Maxwell’s  wave  equation  in  the  Coulomb  gauge, 


I'  1  ^ 

V  c'  df  j 


—  4n  — ' 
A  =  -  —  J, 


(28) 


where  A  is  the  optical  wave  vector  potential  and  is  the  transverse  current  density. 
For  a  helical  undulator,  the  light  emitted  will  be  circularly  polarized,  and  the  wave  vector 
corresponding  to  the  field  (11)  can  be  written  in  rectangular  coordinates  as 

A(z,t)  =  ^^^’^^(sin\t/,cos\|/,0)  (29) 

k 
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where  \\f  =kz-(X)t  +  ^(z,t) .  We  assume  the  field  amplitude  and  phase  vary  slowly  in 
time  over  an  optical  period  [dEj dt  <Kco(t))  and  in  space  over  an  optical 

wavelength  {^dEjdz  <K  kE,  d^  jdz  <K  A:(t) ) .  Then  the  left-hand  side  of  the  wave  equation 
can  be  written 


r  2  1  5'  ^ 

_ , 

1 

"dE 

1 

dE'] 

A 

«2 

+  - 

— 

^  dt^  j 

c 

dt  y 

-2E 

-e- 

1 

+  - 

fcj)  ^ 

dz 

c 

dt  j 

(cos\|/,-sin\|/,0) 


(sin\|/ ,  cos\|/ ,  0)  . 


(30) 


Define  the  unit  vectors  8j  =  (cos\|/ ,-sin\|/ ,0)  and  =  (sin\|/ ,  cos\|/ ,  0) ,  then  project  the 
wave  equation  onto  each  of  these  vectors  to  obtain 

2n 
c 


dz  c  dt 


^  1  fcj) 

ydz  c  dt ) 


(31) 


Now  write  the  total  beam  current  as  a  sum  of  single-particle  currents: 


(32) 


where  is  the  position  of  an  electron  with  transverse  velocity  given  by  (13).  The 
projection  of  onto  Sj  and  83  is 


J ~ 


ecKp  /  ^k^z)) 

Y 

^^^^(sin(\|/  +kf^z)') 


(33) 


where  we  have  assumed  »  1 ,  and  replaced  the  sum  over  all  the  electrons  by  an  average 
(•••)  over  sample  electrons,  weighted  by  the  current  density  p.  Substituting  (33)  into 


(31)  and  using  +(t)  =\j/  +  A:oZ  ,  the  wave  equation  can  be  written  in  complex  form  as 


y  dz  c  dt  j 


Ee^  =  - 


IneKp 
Y 


N- 


(34) 
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If  we  assume  a  long  electron  pulse,  we  can  remove  the  z-dependence  from  this  equation, 
and  it  simplifies  to 

“  =  ^  =  (35) 

where  a  =  |a|e'''’  is  the  dimensionless  optical  field  and  j  is  the  dimensionless  current 
density.  This  is  known  as  the  PEL  wave  equation.  It  has  the  same  form  for  a  linear 
undulator,  if  we  replace  K  in  the  definitions  of  a  and  j  by  K  (  Jq(^)  -  ,  as  before. 

The  FEE  wave  equation  (35)  clearly  demonstrates  that  the  electron  beam  needs  to 
be  bunched  so  that  'm  order  to  achieve  gain.  If  the  electrons  are  bunched  near 
^  « TT  ,  the  field  amplitude  |  a  \  will  increase;  if  the  electrons  are  bunched  near  «  0 ,  the 
amplitude  will  decrease.  Bunching  near  ^  «±7i:/2  will  drive  the  optical  phase  (j),  with 
little  change  in  the  amplitude.  The  wave  equation  also  shows  that  the  gain  depends  on 
the  dimensionless  current  density  y;  for  y  <  tt  there  is  low  gain,  and  for  y  »  tt  there  is 
high  gain. 


5,  Diffraction 

For  a  low  gain  FEE,  the  fundamental  optical  mode  inside  the  laser  cavity  is 
mainly  determined  by  the  cavity  mirrors,  and  is  typically  Gaussian  in  shape.  The  beam 
waist  radius  is  on  the  order  of  a  few  mm,  and  is  usually  located  in  the  middle  of  the 
undulator.  From  there,  the  beam  will  spread  out  by  diffraction,  doubling  in  area  over  a 
distance  Zq  known  as  the  Raleigh  length,  which  is  determined  by  the  curvature  of  the 

mirrors.  If  Zg  is  too  large,  the  mode  waist  becomes  large,  so  there  is  reduced  overlap 
with  the  narrow  electron  beam  (radius  1  mm),  and  less  gain.  If  Zg  is  too  small,  the 

mode  waist  is  small,  but  the  beam  rapidly  expands  near  the  ends  of  the  undulator,  again 
reducing  overlap  and  gain.  The  value  of  Rayleigh  length  that  theoretically  optimizes 
overlap  and  gain  is  Zg  =  l/ ^/l2  «  0.3E ,  where  L  is  the  undulator  length  [15]. 

To  study  the  evolution  of  transverse  optical  modes,  we  use  the  parabolic  wave 
equation  [16]; 
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(36) 


■-v;+  — 
4  0T 


i[x,y,z,x)  =  -{^je 


{x,y,z-T:,x) 


where  =5^+5^,  and  the  x  and  y  eoordinates  are  normalized  to  the  eharacteristie 

mode  radius,  -^LX/n  .  The  operator  -N\Ia  deseribes  the  beam  diffraction.  The 
electron  evolution  is  determined  by  the  pendulum  equation  (18),  which  still  has  the  same 
form,  but  the  field  amplitude  |a|  now  has  x  andy  dependence. 

The  current  density  j  on  the  right-hand  side  of  (36)  drives  both  the  real  and 
imaginary  parts  of  the  optical  field  a  on  the  left-hand  side.  Thus  a  narrow  electron  beam 
can  distort  the  optical  mode  by  changing  the  amplitude  and/or  the  phase  of  the  central 
portion  of  the  wavefront.  In  an  PEL  oscillator  with  modest  gain,  y  «  10  ,  the  steady-state 
resonator  mode  shape  can  be  altered  by  the  FEE  interaction.  In  an  FEE  amplifier  with 
very  large  gain,  7»10%  a  process  known  as  “optical  guiding”  [17]  can  effectively 
contain  the  light  in  the  vicinity  of  the  electron  beam,  continuously  opposing  the  natural 
diffraction. 


D.  WEAK  OPTICAL  FIELDS 
1,  EEL  Integral  Equation 

In  weak  optical  fields,  |a|  <sC7t:  ,  we  use  the  expansion  (20)  for  ^  and  v,  and  the 
wave  equation  (35)  becomes 


while  the  pendulum  equation  (18)  becomes 


2L 


(38) 


Integrating  the  latter  twice  gives 

C  ^‘^(t)  =  £  +a*  .  (39) 
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Substituting  that  result  into  (37)  and  assuming  the  electron  phases  are  evenly  distributed, 
we  obtain 

a(T )  =  —  [  dx'i  dx" a{x ") .  (40) 

'2^  Jo  Jo  \  /vq 

If  the  electrons  have  a  distribution y(^)  of  initial  phase  velocities  v,  =Vo  +  q  about  Vg ,  the 

/•+00 

Fourier  transform  of  this  distribution  is  F(t)  =  e“^  f{q)dq  .  Using  this  “characteristic 

J  -00 

function”  and  the  substitution  r '  =  t  -t  "  we  can  rewrite  (40)  in  the  form 

a(T)  = dx' a(x  -x')  .  (41) 

Switching  the  order  of  integration  and  evaluating  the  inner  integral, 

a(x )  =  -  't  'F(x  a(x  -x ')  .  (42) 

This  is  known  as  the  FEL  integral  equation  [18].  It  is  valid  only  for  weak  optical  fields, 
but  for  both  the  high-gain  and  low-gain  regimes. 

2.  High  Gain 

We’ve  previously  derived  a  formula  for  low  gain  (y  <  tt  )  in  weak  fields  (26).  In 
the  high  gain  regime  (y  »  tt  ) ,  the  optical  field  amplitude  and  phase  evolve  rapidly.  For  a 
perfect  electron  beam,  F{x  )  =  1 ,  (42)  can  be  integrated  to  describe  high  gain  at  resonance 
(Vg  =  0)  as  a  function  of  time  [1], 

G{x)^^exp[43{j/2fx).  (43) 

The  growth  of  the  field  is  exponential  and  is  maximum  at  resonance  for  large y. 

3,  Electron  Beam  Quality  Effects 

The  FEL  integral  equation  (42)  is  useful  for  studying  the  effects  of  electron  beam 
quality  (energy  spread  and  emittance).  Note  that  it  contains  no  reference  to  individual 
electrons,  so  it  can  be  used  to  determine  the  field  evolution  for  a  given  distribution  of 
initial  electron  phase  velocities,  without  having  to  follow  the  motion  of  each  electron  via 
the  pendulum  equation.  Looking  at  the  integrand,  we  see  that  the  term  a(T-T') 
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contributes  to  the  growth  of  the  field  a(T).  This  self-feedbaek  leads  to  exponential 
growth  of  the  optieal  field.  However,  poor  beam  quality  will  eause  the  magnitude  of  the 
eharaeteristie  function  F(t)  to  decay  rapidly,  decreasing  the  growth  rate.  The 
eharaeteristie  function  corresponds  to  the  bunehing  of  the  eleetron  beam  as  it  passes 
through  the  undulator.  As  the  magnitude  of  F(t)  deereases  over  time,  so  does  the 
bunehing  of  the  eleetron  beam  that  gives  rise  to  growth  of  the  optieal  field  in  an  PEL. 
The  eharaeteristie  function  provides  a  visual  means  for  observing  how  different  types  of 
eleetron  distributions  will  affeet  bunehing,  gain,  and  optieal  phase  evolution. 

In  the  early  1990’s,  we  studied  the  effeets  of  various  types  of  eleetron 
distributions  on  the  FEE  meehanism.  Current  beam  diagnosties  in  experiments  do  not  tell 
us  the  preeise  shape  of  these  funetions,  so  it  is  important  to  study  a  variety  of  oases  in 
order  to  understand  the  full  range  of  possible  effeets.  The  author  was  deeply  involved  in 
this  research,  analyzing  the  characteristic  functions  produced  by  various  types  of  electron 
distributions,  and  their  effeot  on  FEE  bunehing  and  gain,  and  was  the  lead  author  of  an 
artiole  we  published  on  this  work.  [12] 

As  an  example,  a  Gaussian  energy  spread  of  rms  width  Ay  leads  to  a  Gaussian 
phase  velooity  distribution 

f{q)  =  (44) 

V  271  CT  Q 

where  <3  q  =  An  N jy  is  the  rms  spread  in  phase  velooities  about  Vg  with  q=v  -Vq. 
This  in  turn  yields  a  Gaussian  eharaeteristie  function. 

Figure  6  shows  that  a  narrow  phase  velooity  spread, in  the  upper-left,  yields  a 
broad  oharaoteristio  funotion,  A  (t  )  in  the  upper-right,  so  bunehing  is  preserved  and  the 
optieal  field  grows  rapidly.  In  oontrast,  a  broad  phase  velooity  spread,  in  the  lower- 
left,  produces  a  narrow  oharaoteristio  funotion,  A  (t  )  in  the  lower-right,  so  bunehing  is 
quiokly  destroyed  and  there  is  muoh  less  gain. 
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Figure  6.  Phase  velocity  distribution  f(q)  and  characteristic  function  F(t)  for 
good  beam  quality  (top)  and  poor  beam  quality  (bottom).  With 
poor  beam  quality,  F(t)  falls  off  faster,  implying  less  bunching  and 
reduced  gain.  (From  [12].) 


To  study  emittance  effects,  we  use  a  distribution  function  p{^x,y,Q where 
(x,y)  is  the  position  of  an  electron  when  it  enters  the  undulator  and  (0^,0 is  the 

injection  angle.  Emittance  effects  would  normally  require  a  3D  simulation  in  order  to 
include  the  changes  in  the  FEE  mechanism  as  the  electrons  drift  off  axis.  However,  it  is 
possible  to  convert  the  distribution  function  to  the  corresponding  phase  velocity 
distribution,  and  then  obtain  the  characteristic  function.  Eor  example,  an  electron  beam 
with  a  linear  spread  has  a  small  angular  spread  in  the  x  direction  only.  The  resulting 
phase  velocity  distribution  is  given  by 

m=  I  (45) 

for  ^<0  and  0  otherwise,  where  ag  =  4nNj^  l{\  + is  the  spread  in  phase 

velocities  induced  by  the  angular  spread  of  rms  width  A0  .  The  amplitude  and  phase  of 
the  characteristic  function  for  this  and  several  other  distributions  is  shown  in  Eigures  7 
and  8. 
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Figure  7.  Amplitude  of  the  eharaeteristie  funetion  F(gt)  for  various  eleetron 
distributions,  as  described  in  the  text.  (From  [12].) 
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Figure  8.  Phase  of  the  characteristic  function  F(aT)  for  various  electron 
distributions,  as  described  in  the  text.  (From  [12].) 
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The  circular  distribution  arises  when  there  is  an  equal  distribution  in  angles  in  the 
X  and  y  directions.  The  triangular  distributions  occur  when  an  emittance  filter  is  used  to 
reject  all  electrons  outside  a  given  radius,  and  outside  a  given  angle;  the  negative  slope 
occurs  with  betatron  focusing,  whereas  the  positive  slope  occurs  if  there  is  no  focusing. 
Notice  in  Figures  7  and  8  that  the  various  distribution  functions  produce  remarkably 
different  characteristic  functions,  and  hence  affect  the  FEL  gain  and  phase  evolution  in 
different  ways.  This  indicates  that  experimentalists  should  pay  attention  not  only  to  the 
width,  but  also  the  shape  of  the  electron  beam  distributions. 


E,  STRONG  OPTICAL  FIELDS 
1.  Saturation 

The  dimensionless  optical  field  strength  |a|  drives  the  electron  bunching 

according  to  the  FEL  pendulum  equation  (18).  Electron  bunching  causes  the  optical  field 
to  grow  via  the  EEL  wave  equation  (35).  As  the  field  grows,  the  peak-to-peak  height  of 

the  separatrix  4|a|  increases,  trapping  a  large  fraction  of  the  electrons  in  the  closed  orbit 

region  of  phase  space,  as  shown  in  Eigure  9.  This  simulation  uses  a  more  realistic 
electron  beam  than  before,  with  1000  sample  electrons  injected  with  a  Gaussian  energy 
spread  =  1  about  the  initial  phase  velocity  Vq  =  5  .  The  initial  electrons  are  shown  in 

yellow.  The  initial  optical  field  strength  is  ag  =  20 .  The  pendulum  and  wave  equations 

are  still  valid  in  strong  fields,  and  the  trapped  electrons  continue  to  follow  the  pendulum 
phase  space  paths.  Halfway  through  the  undulator,  the  electrons  shown  in  light  brown 
are  clearly  bunched  near  the  optimum  phase  for  gain,  ^  =  tt  ,  so  they  give  up  energy  to 
the  optical  wave,  causing  some  electrons  to  lag  behind  others  in  the  beam.  When  they 
reach  the  end  of  the  undulator,  the  electrons  (now  shown  in  red)  are  bunched  at  the 
optimum  phase  for  absorption,  ^  =  0  •  This  “over-bunching”  leads  to  reduced  gain  and 
the  onset  of  saturation. 
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Figure  9.  Electron  evolution  in  (C  ,v  )  phase  space,  optical  field  gain  G(t) 
and  phase  (t  )  for  a  single  pass  through  a  high  power  EEL 
(flo  =  20) .  The  electrons  “over-bunch”  and  the  gain  is  reduced, 
indicating  the  onset  of  saturation. 


At  the  onset  of  saturation,  the  mean  phase  velocity  of  the  electron  bunch  is 
reduced  on  the  order  of  the  peak-to-peak  height  of  the  separatrix,  Av  =  4|a|  ~2tz  .  This 
is  related  to  the  energy  extracted  from  the  bunch  by  Av  =4nN  Ay  jy  ~2n  ,  which  gives  a 
maximum  theoretical  efficiency  r]  s  Ay /y  « 1/2A .  So  a  typical  EEE  with  A  =  100 
periods  has  an  extraction  efficiency  ofri~0.5%.  for  a  high-gain  EEE  {j  »  1),  the  onset 

of  saturation  occurs  at  a  larger  field  strength,  |a| » 2(7'/2)^^^,  and  the  theoretical 
efficiency  is  p  =(^j  12^^  jSN  .  [1] 

2,  Tapering 

To  enhance  the  extraction  efficiency,  it  is  possible  to  taper  either  the  undulator 
period  or  field  strength  along  the  axis,  thus  maintaining  resonance  as  the  electrons  lose 
energy  to  the  optical  field.  This  gives  an  artificial  acceleration  to  the  electron  phase, 
equivalent  to  adding  a  constant  torque  to  the  pendulum  equation,  which  now  becomes 


28 


(46) 


dX 


=  5  +  a  cos(^  +(})) 


where  the  phase  aceeleration  is  5  «  -InN  IsXq/Xq  if  the  undulator  period  is  deereased,  or 


5  «- 


AnN  l{\  +  (AAi/if)  if  the  undulator  field  strength  is  decreased  [1]. 


If  the  taper  rate  is  too  large,  5  >  |a|,  there  will  be  no  electrons  trapped  in  closed 
orbits  near  resonance,  and  there  will  not  be  significant  gain.  If  the  taper  rate  is  too  small, 

I  ll/2 

5  <  4|a|  ,  it  will  not  overcome  the  natural  deceleration  of  the  electrons,  and  there  is  no 

advantage  to  tapering.  Also,  tapering  is  only  useful  in  strong  optical  fields,  |a|>7i:  , 

when  saturation  begins  to  reduce  the  gain.  All  of  these  criteria  for  tapering  can  be 
summarized  in  the  following  expression: 

271  <4|a|‘'^<5<|a|.  (47) 

Tapering  postpones  saturation  and  enhances  the  natural  efficiency  of  an  PEL; 
however,  it  also  reduces  the  weak-field  gain,  so  the  technique  is  more  commonly  used 
with  FEE  amplifiers  rather  than  oscillators.  Negative  tapering  has  also  been  investigated, 
and  has  been  shown  to  be  useful  in  some  cases  [19]. 


3.  Trapped-Particle  Instability 

In  strong  optical  fields,  many  electrons  are  trapped  in  the  closed-orbit  region  of 
phase  space,  as  shown  in  Figure  9.  The  trapped  electrons  begin  to  oscillate  inside  the 
separatrix.  A  perturbation  expansion  of  the  pendulum  equation  near  the  stable  fixed  point 
,v )  =  (71/2,0)  gives  an  equation  for  these  synchrotron  oscillations  [1], 

C  (^)*Co+— sin(vj)  (48) 

V 

s 

1  /7 

where  (^o,Vo)  is  the  initial  position,  and  =  \a\  is  the  synchrotron  frequency.  This 

modulation  of  the  electron  beam  causes  sidebands  to  appear  around  the  fundamental  at 
Vg  ±Vj,  corresponding  to  a  wavelength  shift  AX/X  =vj2nN,  which  is  the  ratio  of  the 
number  of  synchrotron  oscillations  to  the  number  of  undulator  periods.  An  optical  field 
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strength  of  |a|  =  47t:  ^  is  required  for  one  eomplete  oscillation.  High  power  PEL  amplifiers 

can  have  many  tens  of  synchrotron  oscillations,  causing  significant  sideband  growth  and 
fluctuations  of  the  output  power.  In  an  FEE  oscillator,  repeated  low  gain  with  about  one 
synchrotron  oscillation  each  pass  is  sufficient  to  generate  sidebands. 

This  “trapped-particle  instability”  can  destroy  FEE  coherence,  but  also  enhances 
the  total  output  power,  so  it  can  be  desirable  in  some  applications  where  high  power  is 
more  important  than  a  narrow  spectrum.  Sidebands  can  be  suppressed  by  decreasing  the 
current  or  increasing  the  cavity  losses,  or  by  tapering  the  undulator,  which  decreases  the 
fraction  of  trapped  electrons.  In  the  latter  case,  the  synchrotron  frequency  is  decreased  to 

/,  |2  ,\V4 

V.»(H  -8")  ■[!] 

F.  SHORT  PULSES 

RE  linacs  are  capable  of  producing  very  short  electron  pulses,  on  the  order  of 
picoseconds,  which  can  create  very  short  optical  pulses  in  an  FEE.  The  relativistic 
electrons  travel  slightly  slower  than  the  photons.  The  resonance  condition  implies  that 
during  one  pass  through  the  undulator,  the  electron  pulse  will  lag  behind  the  optical  pulse 
by  a  slippage  distance  NX  .  If  the  electron  pulse  length  is  on  the  order  of  or  less  than  the 
slippage  distance,  then  it  will  not  overlap  the  optical  pulse  over  the  entire  pass  through 
the  undulator,  and  gain  will  be  reduced.  The  electrons  at  the  leading  edge  of  the  pulse 
see  the  full  optical  field  strength  and  bunch  accordingly.  As  the  bunched  electrons  slip 
back  relative  to  the  optical  pulse,  they  amplify  the  trailing  edge  of  the  optical  pulse.  Thus 
the  optical  pulse  centroid  travels  slower  than  the  speed  of  light,  a  process  called  “pulse 
lethargy”  [1]. 

1.  Desynchronism 

If  the  electron  pulses  are  timed  so  the  next  pulse  arrives  at  the  beginning  of  the 
undulator  at  the  same  time  as  the  rebounding  optical  pulse  (exact  synchronism),  the 
optical  pulse  will  gradually  fall  behind  the  electron  pulse,  and  optical  power  will  decay. 
To  compensate,  the  optical  pulse  must  arrive  slightly  ahead  of  the  electron  pulse.  In 
practice,  this  is  done  by  adjusting  the  cavity  length  S  until  the  maximum  output  signal  is 


30 


obtained  (“tuning”  the  cavity).  The  artificial  advancement  or  “desynchronism”  of  the 
light  pulse,  relative  to  the  slippage  distance,  is  given  hy  d  =  -2  AS  I  NX  . 

For  large  values  of  desynchronism,  J«0.1,  the  optical  pulse  advances 
significantly  ahead  of  the  electron  pulse  each  pass,  and  only  the  trailing  edge  of  the 
optical  pulse  is  amplified.  This  results  in  a  broad,  low -power  optical  pulse,  with  a  narrow 
spectrum.  The  power  is  generally  too  low  to  trigger  the  trapped-particle  instability.  This 
regime  could  be  desirable  when  stable,  low-power  operation  is  required.  If  d  is  too  large, 
J  »  0.1,  the  electron  and  optical  pulses  won’t  overlap  much,  so  that  gain  will  be  below 
threshold,  and  there  will  be  no  output  power. 

At  moderate  values  of  desynchronism,  J«0.01,  the  optical  field  strength 
becomes  large  enough  to  induce  synchrotron  oscillations,  creating  modulations  of  the 
optical  pulse.  The  desynchronism  causes  these  modulations  to  continuously  drift  from 
the  trailing  to  the  leading  edge  of  the  pulse.  The  optical  power,  spectrum,  and  pulse 
shape  oscillate  periodically  over  hundreds  of  passes;  this  is  known  as  limit  cycle  behavior 
[20]. 

At  smaller  values  of  desynchronism,  J«10\  the  total  output  power  can  be 
larger,  but  the  optical  pulse  can  become  chaotic  due  to  the  trapped-particle  instability.  In 
this  case  the  optical  pulse  is  short  and  modulated  by  sharp  spikes.  [21] 


31 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


32 


III.  SINGLE-MODE  SIMULATIONS 


Prior  to  developing  a  large  comprehensive  simulation  that  covers  all  aspects  of  the 
PEL  interaction,  we  have  for  many  years  used  a  set  of  smaller  simulations  that  model 
various  physical  effects  separately,  enabling  us  to  analyze  the  origins  and  consequences 
of  those  effects.  These  simulations  are  well  understood,  easy  to  modify,  and  can  run  on 
personal  computers.  They  have  been  benchmarked  against  FEE  theoretical  formulas,  and 
by  comparison  to  numerous  FEE  experiments. 

This  chapter  describes  our  single-mode  simulations,  which  model  the  interaction 
of  a  single  optical  mode  (a  single-frequency  plane  wave)  with  the  relativistic  electron 
beam  as  it  travels  through  the  undulator.  Although  higher-order  modes  are  present  in  any 
real  experiment,  we  learn  much  about  the  basic  FEE  physics  from  these  single-mode 
simulations. 

A,  PHASE  SPACE  EVOLUTION 

Our  simplest  simulation,  called  psall,  follows  sample  electrons  as  they  evolve  in 
phase  space  ,v )  and  exchange  energy  with  the  optical  field  over  one  pass  through  the 

undulator.  The  optical  field  can  start  with  a  given  amplitude  and  phase,  or  it  can  be 
allowed  to  develop  from  spontaneous  emission.  The  sample  electrons  are  initially 
distributed  uniformly  in  phase  from  C,  =0  to  2n  ,  but  may  have  a  random  spread  in  phase 
velocities  determined  by  beam  emittance  and  energy  spread.  The  undulator  can  have  a 
positive  or  negative  taper.  In  each  numerical  time  step,  the  pendulum  equation  (18) 
determines  the  microscopic  motion  of  the  electrons  in  phase  space,  while  the  wave 
equation  (35)  is  used  to  self-consistently  update  the  optical  field  amplitude  and  phase. 
The  output  of  psall,  as  shown  in  Figure  9,  displays  the  evolution  of  the  electrons  in  phase 
space,  and  plots  the  optical  field  gain  and  phase  versus  time.  This  simulation  is  useful  for 
understanding  basic  FEE  physical  processes,  such  as  electron  bunching,  optical  field 
growth  and  saturation,  as  discussed  in  the  previous  chapter. 


33 


B. 


GAIN  SURFACES 


We  wrote  another  computer  program  called  gain  which  runs  the  same  single¬ 
mode  simulation  numerous  times,  varying  the  initial  electron  phase  velocity  Vg,  to 
produce  a  gain  spectrum  like  that  shown  in  Figure  4.  In  another  program  called  Gav,  the 
initial  optical  field  strength  is  also  varied  to  produce  a  gain  surface  as  shown  in  Figure 

10.  The  single-pass  gain  is  plotted  versus  initial  electron  phase  velocity  Vg  and  optical 
field  strength  a^,  for  a  low-gain  FEL  (7=1).  The  white  lines  on  the  surface  represent 
contours  of  equal  gain.  Notice  that  as  the  optical  field  strength  is  increased  beyond 
saturation  (ag  >7r),  the  gain  decreases.  The  gain  spectrum  becomes  broader  since  the 

height  of  the  separatrix  increases  with  The  value  of  Vg  that  maximizes  the  gain 

also  increases,  which  corresponds  to  a  longer  wavelength  of  laser  operation;  this  effect 
has  been  observed  in  FEL  experiments  [22]. 
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Eigure  10.  Single-pass  gain  versus  initial  electron  phase  velocity  Vg  and 
optical  field  strength  ,  for  a  low-gain  EEL  (7  =  1).  As 

increases,  the  gain  spectrum  becomes  broader,  the  peak  gain 
decreases  and  shifts  to  a  larger  value  of  Vg . 
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For  a  high-gain  FEL  (7  »  l) ,  the  shape  of  the  gain  surface  changes  significantly. 
Figure  1 1  gives  an  example  of  a  gain  surface  for  an  FEL  with  7  =  100  .  On  the  vertical 
axis  we  plot  the  logarithm  of  the  gain,  ln(l  +  G),  because  the  gain  can  now  change  over 
several  orders  of  magnitude.  Observe  that  in  weak  fields  («(,  <  71 ) ,  the  spectrum  is  now 
broader  (compared  to  Eigure  10)  and  more  symmetric  around  resonance,  Vq  =  0 .  As  the 
field  strength  increases,  the  peak  again  moves  away  from  resonance.  Also  notice  the 
sharp  dip  in  the  surface  near  Vo~-10  in  weak  fields;  since  the  vertical  axis  is 
logarithmic,  this  corresponds  to  gain  «  0  . 
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Eigure  11.  Single-pass  gain,  ln(l  +  G),  versus  initial  phase  velocity  Vg  and 
optical  field  strength  a^,  for  a  high-gain  EEL  (7  =100) .  In  weak 
fields,  the  gain  spectrum  is  broad  and  nearly  symmetric  about  the 
peak  gain  at  Vg  «  0  .  As  ag  increases,  the  peak  gain  decreases  and 
moves  away  from  resonance. 
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1,  Gain  Degradation  Due  to  Electron  Beam  Quality 

Instead  of  varying  the  initial  optical  field,  we  can  vary  other  parameters  such  as 
the  electron  beam  energy  spread  or  an  angular  spread  due  to  beam  emittance.  For 
example,  Figure  12  shows  gain  surfaces  for  circular  and  linear  spreads.  Recall  from  our 
earlier  discussion  that  a  linear  spread  is  an  angular  spread  in  the  x  direction  only,  whereas 
a  circular  spread  in  an  angular  spread  in  both  the  x  andy  directions.  In  Figures  12  (a)  and 
(b)  we  see  that  as  the  angular  spread  Qg  increases,  the  gain  naturally  falls  off,  but  the 
spectrum  details  are  different.  Note  the  different  scales  on  the  Qg  axes.  It  is  clear  that  a 
linear  angular  spread  is  less  susceptible  to  gain  degradation  than  a  circular  spread. 


ln(l+G) 

4^ 


Figure  12(a).  Single-pass  gain,  ln(l  +  G) ,  versus  initial  phase  velocity  Vg 
and  angular  spread  Gg ,  for  a  circular  spread.  As  Gg 
increases,  the  gain  falls  off  rapidly.  (From  [23].) 


36 


Ind+G) 


Figure  12(b).  Single-pass  gain  ln(l  +  G)  versus  initial  phase  velocity  Vg 

and  angular  spread  ag^,  for  a  linear  spread.  As  Gg^ 

increases,  the  gain  falls  off,  but  not  as  rapidly  as  it  does  for 
the  circular  spread.  (From  [23].) 


2,  FEL  Klystron 

Single-mode  simulations  are  also  useful  for  examining  basic  physical  effects  in 
new  undulator  designs.  For  example,  an  FEL  klystron  consists  of  two  identical 
undulators  separated  by  dispersive  magnets  equivalent  to  a  drift  section  of  length 

undulator  periods  [24].  In  the  dispersion  section  between  the  undulators  (at  t  =0.5), 
each  electron’s  phase  evolves  by  =v D ,  where  v  is  the  electron’s  phase  velocity  and 
D  =  NjjN  is  the  dimensionless  dispersive  length.  This  enhances  electron  bunching  and 
thus  the  gain  in  weak  optical  fields.  However  in  strong  optical  fields,  it  precipitates 
overbunching  and  thus  reduces  efficiency.  Figure  13  shows  a  gain  surface  for  a  typical 
FEL  klystron.  Compared  to  the  standard  EEL  shown  in  Eigure  10  for  the  same  current 
density  (7=1),  the  peak  gain  in  weak  fields  is  now  about  three  times  greater,  the  gain  in 

strong  fields  is  significantly  reduced,  and  there  is  more  structure  in  the  gain  spectrum, 
with  multiple  peaks  of  nearly  equal  gain. 
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Figure  13.  Single-pass  gain,  ln(l  +  G),  versus  initial  phase  velocity  Vg  and 


optical  field  strength  a^,  for  an  FEL  klystron  (Z)  =  1) .  Compared 

to  the  standard  FEL  in  Eigure  10,  there  is  now  more  gain  in  weak 
fields,  but  less  gain  in  strong  fields,  and  the  spectrum  has  more 
structure. 


In  the  mid- 1 990 ’s  we  studied  the  SELENE  EEL,  a  ground-based  laser  designed 
for  satellite  power  beaming  [25,26].  This  experiment  used  a  three-stage  klystron 
oscillator  with  dispersive  magnets  at  x  =1/3,  x  =2/3,  and  r  =1,  followed  by  a  single 
pass  EEL  amplifier  referred  to  as  the  radiator.  Our  simulations  showed  the  klystron 
oscillator  strongly  bunches  the  electrons  so  they  begin  spontaneous  radiation  early  in  the 
radiator  and  are  amplified  along  the  entire  length,  with  final  amplified  power  from  the 
radiator  about  500  times  the  level  in  the  oscillator.  The  klystron  configuration  also  limits 
the  saturated  power  within  the  oscillator.  This  scheme  could  produce  a  powerful  output 
beam  without  damaging  the  oscillator  mirrors. 

3,  Positive  and  Negative  Undulator  Tapers 

Last  year  we  published  a  paper  in  which  we  used  single-mode  simulations  to 
compare  positive  and  negative  undulator  tapers  [27].  As  discussed  in  the  above  section 
on  strong  optical  fields,  the  conventional  method  for  enhancing  EEL  efficiency  is  to 
decrease  the  magnetic  field  strength  along  the  undulator  to  maintain  resonance  and 
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extract  more  energy  from  the  electrons.  A  decreasing  magnetic  field  strength, 
AKjK  <0,  adds  a  positive  phase  acceleration 


5  =- 


f  4nNK^^ 

f  AK'] 

[l  +  K^  ) 

1  ^  J 

>0 


(49) 


to  the  pendulum  equation  (46),  so  we  refer  to  this  as  a  positive  taper.  Figure  14  shows  a 
typical  electron  phase  space  after  one  pass  through  an  undulator  with  a  positive  taper, 
5  =  Stt  ,  in  strong  optical  fields,  Oq  =  40 .  This  field  strength  would  normally  be  beyond 
saturation,  but  the  gain  plotted  on  the  right  shows  that  the  field  is  still  growing. 


***  pel  Phase  Space  Evolution  *** 
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Figure  14.  Electron  phase  space,  optical  field  gain  and  phase  evolution  for  an 
FEE  with  a  positive  undulator  taper,  5  =  Sti  ,  in  strong  optical 
fields,  Qq  =  40 .  About  half  of  the  electrons  are  trapped  inside  the 

separatrix,  and  there  is  still  positive  gain,  even  though  the  field 
strength  is  beyond  the  normal  saturation  level. 
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Notice  in  Figure  14  that  the  separatrix  is  shrunk  and  distorted.  About  half  of  the 
eleetrons  are  trapped  in  elosed  orbits  inside  the  separatrix;  the  other  half  are  in  open 
orbits  and  do  not  interact  much  with  the  optical  field.  The  smaller  separatrix  confines  the 
trapped  electrons  near  the  optimum  phase  for  gain,  .  The  fiuetuations  in  the  gain 
are  due  to  synehrotron  oscillations  of  the  trapped  eleetron  bunch  inside  the  separatrix. 

Figure  15  shows  the  phase  spaee  evolution  for  an  undulator  with  a  negative  taper 
of  5  =  -Stt  .  This  corresponds  to  an  inereasing  magnetie  field  strength  along  the 
undulator,  AKjK  <  0 .  The  separatrix  is  now  reversed  from  the  previous  ease,  and  there 
are  no  trapped  eleetrons  inside  it.  However,  the  phase  spaee  paths  of  the  eleetrons  eause 
them  to  buneh  near  =n  about  halfway  down  the  undulator,  as  shown  by  the  light 
brown  dots,  and  they  still  produee  gain. 


***  pel  Phase  Space  Evolution  *** 
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Figure  15.  Eleetron  phase  spaee,  optieal  field  gain  and  phase  evolution  for  an 
FEE  with  a  negative  undulator  taper,  5  =  -Stt  ,  in  strong  optieal 
fields,  Qg  =  40 .  No  eleetrons  are  trapped  inside  the  separatrix,  but 
they  still  buneh  and  produee  gain. 
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Figure  16  shows  FEL  gain  surfaces  for  positive  and  negative  tapers.  The  two 
surfaces  are  related  by  G(vg,d  ,ag)  ~  -G ( -v g , -5 ,  Aq ) .  In  both  cases,  the  weak-field  gain 

is  reduced  but  the  strong-field  gain  is  enhanced  compared  to  the  untapered  case  (Figure 
10).  The  large  tapers  produce  complicated  gain  spectra,  with  multiple  peaks  shifted  away 
from  resonance.  For  the  positive  taper.  Figure  16(a)  shows  that  the  peak  gain  drifts  to 
larger  phase  velocity  Vg,  corresponding  to  a  wavelength  shift  Ak/A,  w  4% ,  as  the  FEL 
evolves  from  weak  fields  to  strong  fields,  for  the  negative  taper,  figure  16(b)  shows  that 
the  peak  gain  remains  at  about  the  same  Vg  as  ag  increases,  so  there  is  no  wavelength 
shift. 


ln(l+G) 

if 


figure  16(a).  Single-pass  gain  In(l-l-G)  versus  initial  phase  velocity  Vg 
and  optical  field  strength  Ug,  for  an  EEL  with  a  positive 

undulator  taper,  5  =  Stt  .  Compared  to  an  untapered  EEL 
(figure  10),  there  is  less  gain  in  weak  fields,  but  more  gain 
in  strong  fields.  As  Ug  increases,  the  peak  gain  shifts  to  a 

larger  value  of  Vg . 
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Figure  16(b).  Single-pass  gain  ln(l  +  G)  versus  initial  phase  velocity  Vg 
and  optical  field  strength  a^,  for  an  FEL  with  a  negative 

undulator  taper,  5  =  -Sn  .  Compared  to  an  untapered  FEL 
(Eigure  10),  there  is  less  gain  in  weak  fields,  but  more  gain 
in  strong  fields.  As  increases,  the  peak  gain  stays  near 

the  same  value  of  Vg . 
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IV.  LONGITUDINAL  MULTIMODE  SIMULATIONS 


In  an  actual  FEL  experiment  with  a  finite  optieal  pulse  length,  there  are  multiple 
longitudinal  modes  inside  the  optieal  eavity.  We  have  developed  several  numerieal 
simulations  that  follow  the  evolution  of  longitudinal  modes  in  an  FEL.  The  optieal  field 
amplitude  and  phase  are  allowed  to  vary  along  the  undulator  (z)  axis.  The  initial  optieal 
field  ean  be  specified  or  allowed  to  develop  from  spontaneous  emission.  The  pendulum 
equation  (18)  and  the  wave  equation  (35)  are  used  to  self-eonsistently  determine  the 
evolution  of  the  eleetrons  and  the  optieal  field  at  eaeh  longitudinal  position.  We  follow 
the  behavior  of  the  optical  pulse  in  a  eo-propagating  window  moving  at  the  speed  of 
light.  By  taking  a  Fourier  transform  of  the  square  of  the  optieal  pulse  amplitude,  we  ean 
obtain  the  optieal  power  speetrum,  and  thus  observe  mode  eompetition  and  eoherenee 
evolution. 

A,  LONG  PULSES 

Consider  an  eleetron  pulse  of  length  4  entering  the  undulator  with  a  eo- 

propagating  optieal  pulse.  In  a  single  pass  through  the  undulator,  the  slower  moving 
eleetron  pulse  slips  baek  relative  to  the  optieal  pulse  by  the  slippage  distanee,  NX  .  If  the 
eleetron  pulse  is  mueh  longer  than  the  slippage  distanee,  =  4/VA.  »  1 ,  a  long  optieal 
pulse  will  form.  In  that  ease,  we  make  the  optieal  pulse  window  a  small,  integral  number 
of  slippage  distanees  long,  and  we  assume  periodie  boundary  eonditions  at  eaeh  end  of 
the  window.  Sinee  the  eleetron  pulse  is  mueh  longer  than  the  optieal  pulse  window,  we 
take  the  eurrent  density  7(2)  to  be  eonstant  in  that  window. 

1,  Single  Pass 

For  studying  long  pulse  FEL  amplifiers,  we  have  a  simulation  ealled  wrapevl  that 
follows  the  eleetron  and  optieal  pulses  over  a  single  pass  through  the  undulator.  The 
output  of  wrapevl  shows  the  evolution  of  the  optieal  pulse,  the  optieal  power  speetrum, 
the  eleetron  energy  distribution,  the  optieal  power  and  gain  versus  dimensionless  time  %  . 
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a,  SLAC X-ray  FEL 

In  the  X-ray  FEL  proposed  at  SLAC  in  1993  [28],  high-energy  7  GeV 
eleetrons  would  be  used  to  produee  radiation  at  a  wavelength  of  4  nm  with  a  peak  power 
of  28  GW.  Since  there  are  no  low-power  input  sources  at  that  wavelength,  and  no 
mirrors  are  available  to  create  an  optical  cavity,  the  LEL  must  create  light  from 
spontaneous  emission  and  amplify  it  to  high  power  over  a  single  pass  through  the 
undulator,  in  a  process  known  as  self-amplification  of  spontaneous  emission  (SASE) 
[29].  Eigure  17  shows  the  output  from  a  simulation  of  this  EEL.  The  high  peak  current 
of  2500  A,  the  60m  long  undulator  with  N  =  723  periods,  and  other  design  parameters 
give  dimensionless  current  density  j  =  2400  ,  so  this  is  a  high-gain  EEL.  The  electrons 
have  a  moderate  energy  spread  g<j  =  3.6  ,  and  a  small  angular  spread  Gg  =  1 .  The  0.16  ps 
duration  electron  pulse  corresponds  to  dimensionless  pulse  length  g^  « 16  . 

In  the  middle  left  of  Eigure  17  we  show  the  evolution  of  the  optical  field 
amplitude  |a(z,T  )|  over  a  single  pass  from  r  =  0  to  1,  as  an  intensity/contour  plot  inside 

a  window  two  slippage  distances  long  (z  =  -l  to  l),  which  moves  with  the  pulse  at  the 
speed  of  light.  The  optical  field  amplitude  is  represented  by  the  color  scale  at  the  top  of 
the  figure,  with  the  largest  amplitude,  |a|  =  362 ,  shown  as  red  and  zero  field  shown  as 

blue.  Observe  that  the  optical  field  grows  from  spontaneous  emission  to  form  a  broad 
pulse  with  fluctuations  due  to  shot  noise,  as  shown  in  blue  in  the  upper  left.  In  the  lower 
left  is  the  evolution  of  the  optical  gain,  ln(l-l-  G(t  )) .  In  the  lower  center  the  theoretical 

weak-field  gain  spectrum  ln(l-l-  G(v  ))  for  7=2400  is  shown  in  black  for  reference.  In 
the  center  of  the  figure  is  an  intensity/contour  plot  of  the  optical  power  spectrum  P{y  ,t  ) , 

which  shows  coherence  developing  over  a  single  pass.  The  final  optical  spectrum,  shown 
in  magenta  in  the  upper-center,  has  a  linewidth  corresponding  to  Ak/A,  «  0.07%  .  The 
triangular  tick  mark  at  the  top  points  to  the  center  of  the  initial  power  spectrum,  and  the 
rectangular  tick  mark  points  to  the  center  of  the  final  power  spectrum.  In  the  lower  right 
is  the  total  optical  power  E’(t)  versus  time,  showing  saturation  near  the  end  of  the 
undulator.  If  higher  power  is  required,  the  undulator  could  be  tapered  to  delay  the  onset 
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of  saturation,  at  the  cost  of  a  broader  linewidth.  In  the  middle-right,  another 
intensity/contour  plot  shows  the  evolution  of  the  electron  phase  velocity  distribution 
/(v  ,T ) .  As  we  discussed  before,  a  change  in  an  electron’s  phase  velocity  corresponds  to 

a  change  in  its  energy,  so  the  final  broad  phase  velocity  distribution,  shown  in  red  in  the 
upper-right,  indicates  that  significant  energy  has  been  extracted  from  the  electrons. 


******  FEL  Wrap  Evolution  ****** 


j=2400  0^=3. 6  <7g=l 
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Figure  17.  Longitudinal  multimode  simulation  results  for  the  SLAC  X-ray 
FEL.  The  various  plots  are  described  in  the  text.  The  optical 
power  spectrum,  T’(v,t),  shows  coherence  evolving  as  the  laser 
line  narrows  over  a  single  pass  through  the  undulator.  (From  [30].) 

2,  Multiple  Passes 

For  studying  long  pulse  FEL  oscillators,  we  have  another  simulation  called 
wrapevn  that  follows  the  electron  and  optical  pulses  over  many  passes  through  the 
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undulator,  again  using  periodic  boundary  conditions  at  each  end  of  the  optical  pulse 
window.  Mirror  transmission  losses  at  the  end  of  each  pass  are  determined  by  the 
resonator  quality  factor  Q.  The  output  of  the  program  shows  the  evolution  of  the  optical 
pulse,  the  optical  and  electron  spectra,  and  the  optical  power  and  gain  plotted  versus  pass 
number  n. 


a.  Boeing  1  kW PEL  Oscillator 

In  the  early  1990’s,  Boeing  designed  an  PEL  oscillator  to  produce  1  kW 
average  power  at  a  wavelength  of  0.6  pm.  It  was  planned  to  be  a  high-gain  FEE, 
7  =  658,  with  N  =  220  undulator  periods,  and  a  fairly  large  initial  electron  energy 
spread,  =11.  The  proposed  resonator  quality  factor  Q  =  5  corresponded  to  ~20% 

losses  per  pass.  Figure  18  shows  the  results  of  a  wrapevn  simulation  of  this  FEE.  This 
figure  is  similar  to  the  previous  examples,  but  now  the  evolution  is  over  n  =  400  passes, 
instead  of  over  a  single  pass.  The  bottom  left  corner  shows  the  gain,  G(n) ,  at  the  end  of 
each  pass  through  the  undulator.  The  FEE  starts  at  resonance  Vq  =0  with  an  initial 
optical  field  of  ag  =  1 .  As  the  electrons  slip  back  past  the  optical  field  sites  and  execute 
synchrotron  oscillations,  they  continually  modulate  the  optical  wave  envelope. 
Numerous  modes  begin  to  grow,  and  the  power  P{n)  increases  dramatically  as 

sidebands  appear  in  the  power  spectrum  P(\  ,n).  The  final  power  spectrum  E’(v )  shows 
several  prominent  sidebands  up  to  5v  «  200  away  from  the  fundamental,  corresponding 
to  about  a  7%  frequency  spread.  The  evolution  of  the  optical  field  amplitude  |a(z,n)|  is 
chaotic,  due  to  nonlinear  interactions  between  the  modes. 
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******  pel  Wrap  Evolution  ****** 
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Figure  18.  Longitudinal  multimode  simulation  results  for  the  Boeing  1  kW 
FEL  oscillator.  The  optical  field  amplitude  |a(z,n)|  shows  the 
development  of  a  chaotic  field  over  n  =  400  passes.  (From  [31].) 


B,  SHORT  PULSES 

If  the  electron  pulse  length  is  just  a  few  slippage  distances  long  or  less  than 
the  slippage  distance,  then  we  must  allow  the  current  density  to  vary  along  the  undulator 
axis.  For  numerical  convenience,  we  choose  a  parabolic  pulse  shape  given  by 
j{z)  =  jil-2z^/ol^  for  7(z)>0.  The  peak  current  density  j  and  electron  pulse 

length  are  specified.  We  follow  the  evolution  of  the  optical  pulse  in  a  co-propagating 
window  long  enough  to  contain  the  entire  pulse,  so  we  do  not  use  periodic  boundary 
conditions  for  these  simulations. 
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1,  Single  Pass 

For  studying  short  pulse  FEL  amplifiers,  we  have  a  simulation  called  pulsevl  that 
follows  the  electron  and  optical  pulses  over  a  single  pass  through  the  undulator.  The 
output  of  the  program  shows  the  evolution  of  the  optical  and  electron  pulses,  the  optical 
and  electron  spectra,  as  well  as  the  optical  power  or  gain  versus  dimensionless  time  %  . 

a,  Boeing  APLE  Amplifier 

In  the  early  1990’s,  we  studied  pulse  slippage  effects  in  the  Boeing 
Average  Power  Laser  Experiment  (APLE),  an  EEL  that  was  designed  to  produce  100  kW 
of  average  power  by  amplifying  a  series  of  short,  intense  pulses  in  a  single  pass  through  a 
tapered  undulator  [32],  This  author  was  directly  involved  in  planning  this  research, 
developing  and  running  the  simulations,  and  interpreting  their  output,  and  was  the  lead 
author  of  a  paper  we  published  on  the  results  of  this  work  [33].  The  450  A  peak  current 
of  each  micropulse  corresponds  to  jF  =  20000 ,  where  j  is  the  peak  dimensionless  current 
density,  and  the  filling  factor  F  =  0.1  is  the  cross-section  area  of  the  electron  beam 
divided  by  the  cross-section  area  of  the  optical  mode.  The  electron  beam  has  a  Gaussian 
energy  spread  =6  and  an  exponential  angular  spread  Gg  =  7 .  The  electron  pulse 

length  is  twice  the  slippage  distance,  g^  =  2.  The  initial  optical  pulse,  produced  by  a 
separate  EEL  oscillator  [34],  has  dimensionless  amplitude  ag  =18  and  a  pulse  length  of 
g^=6  slippage  distances.  The  20%  tapering  rate  gives  a  phase  acceleration  of 
5  =  20071:. 

Eigure  19  shows  the  output  from  a  simulation  of  the  Boeing  APLE  EEL. 
In  the  lower-left,  the  electron  pulse  7(2)  is  shown  in  blue  at  the  beginning  of  the 

undulator  (t  =  O)  and  in  red  at  the  end  of  the  undulator  (t  =  l) .  Over  a  single  pass 

through  the  undulator,  the  slower-moving  electron  pulse  slips  back  relative  to  the  optical 
pulse  by  the  slippage  distance  NX  .  The  window  width  is  4  slippage  distances.  The 
optical  pulse  \a  (z,t  )|  starts  out  broad,  but  the  final  optical  pulse  in  the  upper  left  is  about 

the  same  width  as  the  electron  pulse,  since  only  the  portion  of  the  optical  pulse  that 
overlaps  the  electron  pulse  is  amplified.  The  corresponding  optical  power  spectrum 
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P(v,t)  becomes  broader  over  time.  The  evolution  of  the  electron  phase  velocity 
distribution  /(v  ,t  )  shows  the  taper  starting  at  t  »  0.3  .  Without  the  taper,  the  optical 
field  would  saturate  at  that  point.  The  strong  taper  allows  the  field  to  continue  to  grow, 
as  shown  by  the  total  power  evolution  T’(t).  Tapering  enhances  the  single-pass 
extraction  efficiency  to  r]  »  5% ,  to  achieve  the  100  kW  goal  averaged  over  many  pulses. 
The  final  electron  phase  velocity  distribution  in  the  upper  right  shows  that  about  half  of 
the  electrons  are  trapped  near  resonance  v  =  0  . 


****  PEL  Pulse  Evolution  **** 
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Figure  19.  Longitudinal  multimode  simulation  results  for  the  Boeing  APLE 
amplifier.  The  electron  phase  velocity  distribution  /(v  ,t  )  and  the 

optical  power  evolution  E’(t)  show  the  effect  of  the  linear 
undulator  taper,  5  =  IOOtt  ,  starting  at  r  ~  0.3  .  (From  [33].) 
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2,  Multiple  Passes 

For  studying  short-pulse  FEL  oscillators,  we  have  a  simulation  called  pulsevn  that 
follows  the  electron  and  optical  pulses  over  many  passes  through  the  undulator.  The 
electron  pulse  length  and  the  desynchronism  determine  the  evolution  of  the  optical  pulse. 

Figure  20  shows  the  output  from  a  pulsevn  simulation  of  an  FEL  with  a  short 
electron  pulse,  a^=l,  at  exact  synchronism,  d  =  Q.  The  lethargy  effect  causes  the 

optical  pulse  |a(z,n)|  to  drift  away  from  the  electron  pulse  /(z-t),  and  the  optical 

power  P{n)  decays  over  n  =  1000  passes.  The  electron  phase  velocity  distribution 
/(v  ,n)  is  narrow,  so  very  little  energy  is  extracted  from  the  electron  beam. 
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Eigure  20.  Longitudinal  multimode  simulation  results  for  a  short  pulse  EEL  at 
exact  synchronism  (J  =  0)  over  n  =  1000  passes.  The  optical 

pulse  |a(z,n)|  drifts  away  from  the  electron  pulse  7(z-t),  and 

the  optical  power  E’(n)  decays. 
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Figure  21  shows  a  simulation  with  the  same  parameters,  except  now  there  is  a 
small  desynchronism  (J  =  0.002) ,  which  means  the  cavity  is  slightly  shortened  to  allow 

the  optical  pulse  to  “keep  up”  with  the  electron  pulse.  The  optical  power  evolution  P{n) 
shows  that  the  FEL  now  reaches  steady-state  saturation  after  about  n  =  800  passes.  In 
steady-state,  the  gain  each  pass  equals  the  loss  each  pass  (1/0  =  2%) .  When  the  optical 

field  saturates,  the  bunched  electrons  begin  to  execute  synchrotron  oscillations  as  they 
slip  back  over  the  optical  field  sites.  This  induces  a  slight  modulation  of  the  final  optical 
pulse  in  the  upper  left,  and  the  optical  power  spectrum  P{y  ,n)  becomes  broader. 
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Figure  21.  Longitudinal  multimode  simulation  results  for  a  short  pulse  FEL 
with  small  desynchronism  ((i  =  0.002),  showing  steady-state 

saturation  of  the  optical  power  F’(n) . 
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Figure  22  shows  another  simulation  with  large  desynchronism  i^d  =  0.02) .  Now 

the  optical  pulse  advances  ahead  of  the  electron  pulse  each  pass,  becoming  very  broad,  as 
shown  by  the  final  optical  pulse  in  the  upper  left.  If  the  window  was  wide  enough  to 
contain  the  entire  pulse,  it  would  show  that  the  pulse  has  a  long  exponential  tail  [1].  Less 
overlap  between  the  electron  and  optical  pulses  yields  less  gain,  and  the  steady-state 
power  P{n)  is  reduced.  The  steady-state  optical  power  is  reduced  at  large 

desynchronism,  and  for  large  enough  d,  the  power  is  zero.  This  defines  a  range  of 
desynchronism  values  where  the  FEL  operates. 
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Figure  22.  Longitudinal  multimode  simulation  results  for  a  short  pulse  FEL 
with  large  desynchronism,  (J  =  0.02).  The  optical  pulse  |a(z,n)| 

is  now  long,  and  the  steady-state  power  P{n)  is  reduced. 
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a. 


Limit  Cycles:  Stanford  FIREFLY  PEL 


In  strong  optical  fields,  the  trapped-particle  instability  can  lead  to 
modulations  in  the  optical  pulse.  These  can  appear  as  sharp  spikes  that  drift  periodically 
across  the  pulse  envelope  at  a  rate  depending  on  the  desynchronism.  This  “limit-cycle” 
behavior  has  been  observed  in  PEL  experiments  [35,36].  Figure  23  shows  the  output 
from  a  simulation  of  the  Stanford  FIREFLY  [37],  an  FEE  that  was  designed  in  1994  to 
produce  1  kW  average  power  at  a  wavelength  of  7.6  pm.  Limit  cycles  are  clearly  seen  in 
the  evolution  of  the  optical  field  |a(z,n)|  and  the  total  optical  power  E’(n) .  Sidebands 

appear  in  the  optical  power  spectrum  P{y  ,n),  and  their  separation  roughly  corresponds 

to  the  synchrotron  frequency,  . 
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Figure  23.  Longitudinal  multimode  simulation  results  for  the  Stanford 
EIREFLY  experiment.  The  trapped-particle  instability  causes 
limit-cycle  behavior  in  the  optical  pulse  \a  (z,n)|  and  the  formation 

of  sidebands  in  the  optical  power  spectrum  /"(v ,  n) .  (From  [37].) 
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b.  Ultra-short  Pulse  FELs 

In  the  mid- 1 990 ’s  we  studied  ultra-short  pulse  PEL  oseillators,  sueh  as 
FELIX  in  the  Netherlands  [38],  which  uses  picosecond  electron  pulses.  This  author 
played  a  key  role  in  planning  this  research,  developing  and  running  the  simulations,  and 
interpreting  their  results,  and  was  the  lead  author  of  a  paper  we  published  on  the  results 
of  this  work  [39].  Figure  24  shows  the  results  of  a  simulation  of  an  ultra-short  pulse  FEE 
with  =  0.1  and  very  small  desynchronism  d  =  0.0002  over  n  =  5000  passes. 


****  fel  Pulse  Evolution  **** 

j=20  d=0.0002  Q=100 

2 


|a(z,n)|  P(v,n)  f(v,n) 


-2  z  2  -101  V  101  0  n  5000 


Figure  24.  Longitudinal  multimode  simulation  results  for  an  ultra-short  pulse 
FEL  =0.1),  showing  the  formation  of  a  narrow  optical  pulse 

|a(z,n)|  and  abroad  optical  power  spectrum  P{y  ,n) .  (From  [39].) 


54 


In  Figure  24,  the  final  optieal  pulse  |a(z,n)|  is  about  the  same  length  as 

the  electron  pulse  7(z-t),  and  much  shorter  than  the  slippage  distance.  Since  the 

optical  pulse  is  so  short,  the  electron  pulse  does  not  “see”  much  of  the  light  until  it  is 
more  than  halfway  down  the  undulator,  leaving  less  time  for  bunching  and  gain. 
However,  the  large  dimensionless  current  density  j  =  20  enables  bunching  and  gain  to 

develop  quickly.  The  1%  loss  each  pass  (g  =  100)  allows  the  optical  field  to  grow  to  a 
large  amplitude  over  n  =  5000  passes. 

c.  Electron  Beam  Energy  Modulation:  Stanford  SCA  EEL 

In  1990,  the  Stanford  Superconducting  Accelerator  (SCA)  FEL  was  used 
to  study  the  effects  of  electron  beam  energy  modulation  [40],  which  is  equivalent  to 
oscillation  of  the  electron  phase  velocity  through  the  relationship  Av  =AtiN  Ay  /y  . 
Figure  25  shows  the  output  from  a  simulation  of  the  Stanford  SCA  FEL  over  n  =  2000 
passes,  with  a  sinusoidal  variation  added  to  the  electron  phase  velocity  on  each  pass, 
V  (n)  =Vo sin(27i:  n/Ag) .  The  amplitude  of  the  phase  velocity  modulation  is 
Ag  =2n  ,  and  the  period  is  Ng  =  500  passes.  The  energy  modulations  can  be  seen  in  the 
evolution  of  the  electron  phase  velocity  distribution  /(v,n),  and  the  optical  power 
spectrum  E’(v,n).  As  strong  optical  fields  develop,  mode  competition  narrows  the 

spectrum,  making  it  more  resistant  to  modulation.  This  prevents  the  peak  optical  power 
from  following  the  sinusoidal  modulations  fully  to  the  right,  so  we  observe  flattening  of 
the  P(y,n)  oscillations  for  v  >  0 .  The  total  optical  power  P{n)  also  oscillates  at  the 
same  frequency;  the  dips  in  the  peaks  correspond  to  narrowing  of  the  peaks  in  the  optical 
pulse  evolution,  |a(z,n)|. 
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Figure  25.  Longitudinal  multimode  simulation  of  the  Stanford  SCA  FEL  with 
modulated  electron  beam  energy.  Oscillations  appear  in  the  optical 

pulse  |a(z,n)|,  the  power  spectrum  P{y,n),  the  phase  velocity 
distribution  f{y,n),  and  the  optical  power  P{n) .  (From  [41].) 


C.  CURRENT  NFS  RESEARCH:  TJNAF  FEL 

For  the  past  several  years,  we  have  been  involved  in  design  studies  for  the 
Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  FEL  in  Virginia  [7].  Their 
accelerator  provides  0.5  ps  duration  electron  pulses  at  34.5  MeV  energy  with  a  peak 
current  of  50  A,  at  a  pulse  repetition  rate  of  18.7  MHz.  These  values  correspond  to 
dimensionless  current  density  y  =10  and  electron  pulse  length  cj^  =1 .  Their  undulator 
has  A  =  41  periods  of  length  Lq  =  2.7  cm  each,  with  undulator  parameter  K  =  0.98  .  The 
EEL  produces  up  to  2  kW  of  average  power  at  a  wavelength  of  A,  =  6  pm  inside  a  g  =  10 
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resonator.  It  is  currently  being  upgraded  to  10  kW  [8]  by  increasing  the  electron  beam 
energy,  the  peak  current,  and  the  pulse  repetition  rate.  A  100  kW  PEL  [42]  is  also  in  the 
design  stage. 

In  the  TJNAF  FEE,  after  the  electron  beam  exits  the  undulator,  it  is  recirculated 
back  to  the  accelerator  for  energy  extraction.  Recirculation  enhances  system  efficiency, 
and  reduces  radiation  shielding  requirements  for  the  beam  dump,  but  it  requires  a  small 
induced  electron  energy  spread  (<10%)  for  beam  transport.  A  high  power  FEE  can 

induce  a  significant  electron  energy  spread.  Various  undulator  designs  have  been 
considered  to  optimize  steady-state  optical  power  while  minimizing  the  exhaust  energy 
spread  of  the  electron  beam. 

1.  Positive  and  Negative  Undulator  Tapers 

East  year  we  published  a  paper  on  positive  and  negative  undulator  tapers  [27], 
which  were  explored  experimentally  using  the  TJNAF  FEE  [43].  We  simulated  linear 
taper  rates  of  AKjK  =  ±5%  ,  ±7.5%,  and  ±10%,  corresponding  to  phase  accelerations  of 
5  =  ±471  ,  ±671  ,  and  ±871  in  the  pendulum  equation  (46).  For  each  value  of  the  undulator 
taper  rate,  we  varied  the  desynchronism  from  J  =  0  to  0.4  to  determine  the  optimum 
weak- field  gain  and  steady-state  power.  For  most  of  the  range  of  d,  a  positive  taper 
provides  no  advantage  over  the  untapered  case.  However,  for  a  small  negative  taper 
5  =  -471:  ,  our  simulations  predict  more  steady-state  power  than  the  untapered  case  up  to 
d  -  0.18  .  Negative  taper  also  induces  slightly  less  energy  spread  on  the  electron  beam. 

For  some  values  of  ii<0.1,  we  observe  limit-cycle  behavior.  Figure  26  shows 
the  results  of  a  pulsevn  simulation  of  the  TJNAF  FEE  over  n  =  2000  passes,  with  a  large 
positive  taper  5  =  871:  ,  and  desynchronism  d  =  0.026  .  To  include  shot  noise,  the  electron 
phases  are  given  an  initial  random  spread  of  amplitude  8C,  =  0.0001 .  In  the  evolution  of 

the  optical  power  spectrum  P{y  ,n),  at  about  n  =  200  passes,  the  FEE  suddenly  shifts  to 

a  different  frequency.  This  corresponds  to  the  onset  of  strong  optical  fields,  where  the 
peak  in  the  gain  curve  shifts  to  a  larger  value  of  v  (refer  to  Figure  16(a)).  After  that, 

limit  cycles  appear  as  oscillations  of  the  optical  field  |a(z,77)|,  the  optical  spectrum 
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P{y  ,n),  and  the  electron  phase  velocity  distribution  f{y,n).  The  total  optical  power 
P{n)  is  modulated  by  about  50%,  at  the  same  frequency.  This  limit-cycle  behavior  is 
also  seen  with  negative  tapers. 


****  pel  Pulse  Evolution  **** 
j=10  0=1  Q=10  S=8z  d=0.026  8^=0.0001 


-2  z  2-50  V  50  0  n  2000 


Figure  26.  Longitudinal  multimode  simulation  results  for  the  TJNAF  FEL 
with  large  positive  taper,  5  =  871  .  Limit-cycle  behavior  is 

observed  in  the  evolution  of  the  optical  pulse  |a(z,n)|  and  the 
optical  power  spectrum  F’(v,n).  The  total  optical  power  F’(n)  is 
modulated  by  -50%.  (From  [27].) 
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2,  Step-tapered  Undulator 

Last  year  we  studied  the  use  of  a  step-tapered  undulator  for  the  proposed  100  kW 
TJNAF  PEL.  This  author  played  a  key  role  in  developing  and  running  the  simulations 
and  interpreting  the  results,  and  was  the  lead  author  of  a  paper  on  this  work,  recently 
accepted  for  publication  [44],  Unlike  the  linear  tapered  undulators  discussed  previously, 
the  step-tapered  undulator  would  have  an  abrupt  change  in  the  magnetic  field  halfway 
through  the  undulator.  In  an  experiment,  this  would  be  done  by  having  two  separate 
undulator  segments,  each  with  a  different  gap  between  the  magnets.  The  modification  to 
the  pendulum  equation  is  similar  to  the  linear  taper  (46): 
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T  -- 
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A  +  acos(^  +(})) 
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where  0  ( z )  =  <  is  the  step  function,  and 
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is  the  phase  acceleration  at  t  =1/2  due  to  the  step  taper  AKjK  . 
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Figure  27  shows  the  results  of  a  series  of  longitudinal  multimode  simulations  of 
the  proposed  100  kW  FEE  for  a  negative  step  taper  A  =  -7t:  ,  with  the  desynchronism 
varied  from  J  =  0  to  0.3  . 
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Figure  27.  Extraction  efficiency  r]  versus  desynchronism  d  for  the  TJNAF 
100  kW  FEE  with  negative  step-tapered  undulator,  A  =  -7t:  .  To 
produce  100  kW,  an  efficiency  of  r\  >0.7%  is  required,  which  is 
achieved  for  0.03  <  d  <  0.07 . 


In  Figure  27,  the  extraction  efficiency  p  =  Ay  /y  =  Av/dii:  A  is  determined  from 

the  shift  in  the  average  electron  phase  velocity,  Av  .  Note  that  p  is  just  the  single-pass 
extraction  efficiency;  the  overall  system  “wall  plug”  efficiency  is  enhanced  by 
recirculation  of  the  electron  beam.  The  steady-state  output  power  is  directly  proportional 
to  p  .  Figure  27  shows  the  general  features  of  desynchronism  (or  “cavity  detuning”) 
curves  that  are  commonly  observed  in  our  simulations  as  well  as  in  FEE  experiments 
[43].  At  exact  synchronism,  d  =  0,  there  is  no  output  power,  as  discussed  previously.  At 
small  desynchronism,  d  «0.01,  there  is  a  sharp  rise  in  the  output  power,  peaking  at 
d  «  0.05  in  this  case.  Fimit-cycle  behavior  is  often  seen  in  this  region,  as  shown  in 
earlier  figures.  As  d  is  further  increased,  the  power  slowly  drops  off,  until  eventually 
there  is  no  output  power,  for  d  >  0.24  in  this  case.  For  stable  operation,  it  may  be 
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desirable  to  operate  the  PEL  at  a  value  of  d  slightly  beyond  the  peak  of  the 
desynchronism  curve.  Otherwise,  small  fluctuations  in  the  cavity  length  could  cause 
large  fluctuations  in  the  output  power,  unless  active  feedback  is  employed  to  keep  the 
cavity  length  constant. 

To  achieve  the  goal  of  100  kW  average  power,  an  efficiency  of  r\  >0.7%  is 
required  for  this  FEE.  Figure  27  shows  that  with  a  step-tapered  undulator  A  = -n  ,  the 
100  kW  goal  is  reached  for  desynchronism  values  0.03  <  d  <  0.07  .  For  no  taper,  the 
peak  output  power  is  slightly  reduced,  and  the  range  of  desynchronism  values  that 
produce  the  100  kW  goal  is  narrower,  0.03  <  J  <0.06.  The  induced  electron  energy 
spread  is  similar  in  both  cases.  Ay  /y  «  4%  ,  well  below  the  requirement  for  recirculation. 
Other  values  of  step-taper  that  we  studied,  A  =  tt  ,  ±2n  ,  did  not  reach  the  100  kW  goal. 

Previously  published  results  based  on  simulations  and  experimental  results  from 
CEIO  and  FEEIX  showed  as  much  as  a  75%  efficiency  enhancement  with  a  negative 
step-tapered  undulator  [45-47].  In  contrast,  we  found  only  a  slight  improvement  in 
efficiency  with  a  small  negative  taper  for  the  parameters  investigated.  The  results  of  our 
research  show  that  step-taper  cannot  always  be  as  effective  as  found  earlier,  and  that  an 
FEE  must  be  far  into  strong-field  saturation  before  tapering  can  extend  the  saturation 
limit. 


3.  FEL  Klystron 

East  year  we  published  a  paper  in  which  we  investigated  the  use  of  an  FEE 
klystron  configuration  for  the  TJNAF  10  kW  FEE  [48].  This  FEE  will  use  an 
electromagnetic  undulator,  so  by  changing  the  magnetic  field  strength,  the  undulator 
parameter  can  be  varied  from  X  =  1  to  4.  Thus  according  to  the  resonance  condition  (14), 
the  FEE  wavelength  can  be  tuned  from  2.3  pm  to  20  pm.  Since  the  undulator  consists  of 
two  separate  sections,  it  is  possible  to  have  a  dispersive  section  in  between.  This  would 
enhance  the  gain  in  weak  fields,  but  reduce  the  saturation  power  and  hence  the  efficiency. 
Our  simulations  showed  that  with  no  klystron  (Z)  =  O)  and  undulator  parameter  X  =  1  or 

2,  the  TJNAF  FEE  could  produce  a  maximum  of  14.5  kW  average  power.  With  X  =  3  or 
4,  the  weak-field  gain  is  below  threshold,  and  there  is  no  output  power,  due  mainly  to  a 
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smaller  filling  factor  F  at  long  wavelengths.  Introducing  a  weak  klystron,  D  =  0.25  , 
reduces  the  average  power  to  slightly  above  10  kW  for  K  -  \  or  2,  but  provides  enough 
gain  to  enable  operation  with  if  =  3  at  an  output  power  of  7.5  kW.  With  a  stronger 
klystron,  Z)  >  0.5  ,  all  four  undulator  parameters  produce  some  output  power,  but  less 
than  the  desired  value  of  10  kW. 
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V.  TRANSVERSE  MULTIMODE  SIMULATIONS 


The  simulations  described  so  far  have  all  assumed  plane  waves,  but  in  an  actual 
PEL  experiment,  an  optical  wavefront  with  finite  transverse  dimensions  will  have 
multiple  transverse  modes.  To  study  these  transverse  modes,  we  have  developed  several 
three-dimensional  simulations  that  include  diffraction  and  self-consistent  evolution  of  the 
optical  wavefront  in  x,  y,  and  t  .  These  simulations  do  not  assume  axial  symmetry;  the 
electron  beam  is  allowed  to  have  separate  emittance  in  the  x  andy  directions.  Transverse 
distances  are  normalized  to  the  characteristic  optical  mode  radius  Wq  =  .  The 

parabolic  wave  equation  (36)  is  used  to  follow  the  evolution  of  the  optical  field. 

A,  SINGLE  PASS 

For  simulating  FEE  amplifiers,  we  have  a  three-dimensional  simulation  called 
3dl  that  follows  multiple  transverse  optical  modes  over  a  single  pass  through  the 
undulator.  The  output  of  3dl  shows  the  evolution  of  the  optical  wavefront  and  the 
driving  current,  the  electron  phase  velocity  distribution,  the  optical  power  and  gain,  and 
the  optical  phase  versus  dimensionless  time  t  . 

1.  Optical  Guiding 

As  mentioned  earlier,  the  FEE  interaction  can  distort  the  optical  mode.  If  the 
current  density  is  very  large,  7»10\  a  phenomenon  known  as  “optical  guiding”  can 
contain  the  light  in  the  vicinity  of  the  electron  beam.  Enhanced  overlap  between  the 
electron  beam  and  the  optical  mode  increases  the  FEE  gain.  There  are  two  processes  that 
contribute  to  optical  guiding.  A  narrow,  high-current  electron  beam  can  continuously 
amplify  the  central  portion  of  the  wavefront,  effectively  narrowing  the  optical  mode  area 
in  a  process  known  as  “gain  guiding”.  A  high-current  beam  can  also  continuously  drive 
the  optical  phase,  opposing  the  phase  shift  due  to  natural  diffraction  and  focusing  the 
light  back  towards  the  electron  beam  in  a  process  known  as  “phase  guiding”.  In  the 
following  example,  both  types  of  optical  guiding  are  observed. 
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a. 


Boeing  ABLE  Amplifier 


Figure  28  shows  the  results  from  a  simulation  of  the  Boeing  APLE 
amplifier,  discussed  earlier.  The  box  in  the  upper  right  shows  the  dimensionless 
parameters  for  this  experiment.  The  current  density  is  7  =  2x10%  with  an  electron  beam 
radius  of  =  0.09  .  The  electrons  have  a  Gaussian  energy  spread  of  width  <3^=6  and 
an  exponential  angular  spread  of  width  Gg  =  7 .  The  initial  electron  phase  velocity 
Vg  =15  is  determined  from  FEE  gain  surfaces  for  these  parameters.  The  initial  optical 
field  amplitude  is  ag  =  1 8  .  The  position  of  the  optical  waist  is  at  the  beginning  of  the 
undulator,  =  0,  in  an  FEE  amplifier  with  no  resonator  mirrors.  The  dimensionless 
Rayleigh  length,  Zg  =  Wg  =  0.8  ,  is  determined  by  the  optical  mode  radius,  Wg .  The 

undulator  has  N  =  256  periods,  with  a  20%  positive  taper  rate  corresponding  to  a  phase 
acceleration  of  5  =  2007t:  .  The  color  scale  shown  in  the  parameter  window  is  used  for 
the  intensity  plots  of  the  optical  field  amplitude  and  the  electron  beam  current  density, 
with  dark  blue  corresponding  to  an  amplitude  or  density  of  zero,  and  light  blue 
corresponding  to  the  maximum  amplitude  or  density.  The  white  contour  line  on  the 
intensity  plots  corresponds  to  half  of  the  maximum  value,  as  indicated  on  the  color  scale. 

The  intensity/contour  plot  in  the  upper-left  of  Figure  28  shows  a  slice 
through  the  optical  wavefront,  |a  (x,t  )|  at  y  =  0 ,  as  it  evolves  in  time,  and  next  to  it  is  the 

final  wavefront  |a(x,y)| .  The  narrow,  high-current  electron  beam  causes  optical  guiding, 

which  yields  a  much  smaller  final  spot  size  than  would  be  predicted  from  natural 
diffraction.  The  slow  evolution  of  the  optical  phase  in  the  center  of  the  wavefront, 
(j)(0,T)  in  the  center-right,  indicates  that  the  predominant  process  in  this  case  is  gain 
guiding.  In  the  center-left  is  shown  the  evolution  of  the  driving  current 
G  (x,t)  =  7^cos(^  +^))  y  =  0,  where  j  is  the  dimensionless  current  density  at  the  site 

(x,  y)  and  (•  •  •)  indicates  an  average  over  all  the  electron  phases  at  that  site. 
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Figure  28.  Three-dimensional  simulation  results  for  the  Boeing  APLE 
amplifier.  The  various  plots  are  deseribed  in  the  text.  The 

evolution  of  the  optieal  wavefront  |fl(x,T)|  shows  optieal  guiding. 

(From  [33].) 


The  final  electron  beam  profile  a  (x,  j)  in  the  center  of  Figure  28  and  the 
final  wavefront  |a(x,  j)|  above  it  are  used  to  estimate  the  filling  factor,  F  =  jw^  »  0.1 , 
where  is  the  electron  beam  radius  and  is  the  mean  optical  mode  radius.  In  the 
lower-left  is  the  evolution  of  the  electron  spectrum  /(v,t),  showing  the  20%  taper 
starting  at  the  saturation  time  »  0.2  .  The  final  phase-space  plot  /(^  ,v )  in  the  lower- 
center  shows  that  about  70%  of  the  electrons  remained  trapped  in  closed  orbits.  The 
graphs  in  the  lower-right  show  that  the  optical  power  F’('r)  and  gain  (j('r)  continue  to 

grow  linearly  after  ,  so  the  taper  is  effectively  increasing  the  efficiency.  This 
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simulation  predicts  an  efficiency  of  r]  =  8%  ,  which  is  too  optimistic  since  pulse  slippage 
effects  are  neglected. 


B,  MULTIPLE  PASSES 

For  simulating  transverse  optical  modes  in  FEL  oscillators,  we  have  another 
three-dimensional  simulation  called  3dn  that  runs  for  many  passes,  including  the  effects 
of  mirror  transmission  and  edge  losses.  The  shape  and  size  of  the  cavity  mirrors  is 
determined  by  the  Rayleigh  length,  the  optical  mode  waist,  and  the  mirror  edge  losses. 
The  mirror  transmission  can  be  varied,  as  well  as  the  output  coupling  mechanism.  The 
output  of  3dn  shows  the  evolution  of  the  optical  wavefront  over  many  passes,  the  final 
optical  mode  shape,  the  electron  phase  velocity  distribution,  the  optical  power  and  gain 
versus  pass  number  n. 

1.  Optical  Mode  Distortion 

In  a  typical  FEL  oscillator,  the  steady-state  optical  mode  has  a  transverse 
Gaussian  profile,  determined  mainly  by  the  resonator  geometry  (the  size,  shape  and 
location  of  the  mirrors,  and  the  outcoupling  mechanism).  However,  a  moderate  electron 
beam  current,  y  «  100  to  1000,  can  distort  the  optical  mode,  affecting  the  average  filling 
factor  and  hence  the  weak-field  gain  and  steady-state  power. 

a.  SELENE  EEL  Klystron  Oscillator 

The  SELENE  experiment,  mentioned  previously,  is  a  two-stage  EEL 
proposed  in  1993  for  satellite  power-beaming  applications  [49],  based  on  the  Novosibirsk 
racetrack  microtron  [50].  The  three-stage  EEL  klystron  oscillator  consists  of  three 
identical  undulators  separated  by  two  dispersive  magnets  located  at  dimensionless  times 
T  =  1/3  and  t  =  2/3 .  The  klystron  design  was  intended  to  enhance  electron  bunching  and 
reduce  the  saturated  optical  power  inside  the  cavity.  The  purpose  of  the  oscillator  is  only 
to  bunch  the  electrons  before  they  are  injected  into  an  external  radiator  [25].  The  radiator 
creates  high  average  power  over  a  single  pass  by  self-amplification  of  spontaneous 
emission  from  the  bunched  electron  beam. 
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Figure  29  shows  the  output  from  a  3dn  simulation  of  the  SELENE 
oseillator.  The  dimensionless  parameters  for  this  experiment  are  shown  in  the  upper- 
right.  All  longitudinal  distances  are  normalized  to  the  undulator  length  L,  and  transverse 

distances  are  normalized  to  the  characteristic  optical  mode  radius  Wf^  =  ^LX/ti  .  The 
micro tron  provides  an  electron  beam  with  current  density  j  =  2500  and  radius 
Qg  =  0.14  .  The  electron  beam  can  be  given  an  offset  from  the  undulator  axis,  but  in  this 
case  there  is  no  offset,  so  Xq  =  0 .  At  each  pass,  the  electrons  are  injected  with  an  initial 
Gaussian  spread  in  phase  velocity  of  width  G(j  =  3 ,  centered  on  Vq  =  1 .  The  simulation  is 
started  with  a  weak  optical  field,  aQ=\,  but  the  optical  power  quickly  saturates  after 
n  =  10  passes.  The  Rayleigh  length  Zq=0.4  and  the  position  of  the  optical  waist 
=  0.5  are  chosen  to  maximize  the  overlap  between  the  electron  beam  and  the  optical 
mode.  The  klystron  dispersive  strength  is  given  by  D  =  0.66 .  The  three  undulator 
sections  have  a  total  of  N  =  120  periods.  The  mirror  radius  is  =1.9,  and  radius  of 

curvature  =1.2.  Cavity  losses  are  given  by  0  =  10  corresponding  to  approximately 
10%  mirror  transmission  and  edge  losses  around  the  mirrors  of  10%  per  pass. 

In  the  upper  left  of  Eigure  29,  we  show  a  slice  through  the  optical 
wavefront,  |a(x,n)|  at  y  =  0,  as  it  evolves  over  n  =  10  passes,  and  next  to  it  the  final 

optical  wavefront  |a(x,y)| .  Below  that  is  a  slice  through  the  optical  wavefront,  |a(x,T  )| 
at  y  =  0 ,  as  it  evolves  over  the  final  pass,  with  the  mirrors  shown  at  each  end,  and  tick 
marks  beneath  indicating  the  beginning  of  the  undulator  (r  =  O)  and  the  end  of  the 
undulator  (t  =  l) .  The  mirror  separation  is  artificially  shortened  to  twice  the  undulator 
length  for  numerical  convenience. 
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Figure  29.  Three-dimensional  simulation  results  for  the  SELENE  FEE 
oseillator,  showing  mode  distortion  in  the  evolution  of  the  optical 

wavefront  |a(x,n)|.  (From  [51].) 


Figure  29  shows  optical  mode  distortion,  and  the  effect  of  the  klystron 
magnets  is  clearly  seen  in  the  final  mode  shape,  |a(x,T)|.  In  the  lower-left  is  the 

evolution  of  the  electron  spectrum  /  (v  ,t  ) .  The  final  electron  phase-space  plot  in  the 

lower-center  shows  evidence  of  bunching.  The  achromatic  bending  magnet  used  to 
extract  the  electrons  from  the  oscillator  [52]  would  add  more  dispersion  and  further 
improve  the  bunching.  In  the  lower  right  is  the  evolution  of  the  optical  power  E’(n)  and 

the  gain  G(n),  showing  that  the  oscillator  quickly  saturates.  It  may  be  necessary  to 
operate  the  resonator  at  a  low  0  <  10  to  prevent  the  power  from  growing  excessively  and 
damaging  the  mirrors. 
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b.  TJNAFlOkWFEL 

Figure  30  shows  a  three-dimensional  simulation  of  the  TJNAF  10  kW 
FEL  oscillator  over  n  =  500  passes.  The  dimensionless  parameters  for  this  experiment 
are  shown  in  the  upper-right.  Their  accelerator  provides  an  electron  beam  with  a  current 
density  of  7  =  165  and  a  radius  of  =0.12.  At  each  pass,  the  electrons  are  injected 
with  an  initial  Gaussian  spread  in  phase  velocity  of  width  g<j  =  0.5  ,  centered  on  Vg  =  9  . 
The  Rayleigh  length  is  Zg  =0.4  with  the  optical  waist  located  in  the  center  of  the 
undulator,  t„=0.5.  The  undulator  has  A  =  24  periods.  The  dimensionless  mirror 
radius  is  =  3.2 ,  and  the  dimensionless  radius  of  curvature  is  =  1.2 .  Cavity  losses  are 
given  by  g  =  6  corresponding  to  approximately  15%  mirror  transmission  and  edge  losses 
around  the  mirrors  of  0.1%  per  pass.  The  simulation  was  started  with  a  weak  optical 
field,  flg  =0.1,  and  allowed  to  develop  until  steady-state  power  was  obtained  in  a  stable 
resonator  mode. 

In  weak  fields,  a<Ti  ,  the  narrow  electron  beam  distorts  the  resonator 
mode,  as  seen  in  the  upper-left  plot  |a(x,n)|  in  Figure  30.  Mode  distortion  reduces  the 

filling  factor,  enhancing  the  gain  in  this  example  by  about  a  factor  of  4  compared  to  a 
simple  theoretical  estimate.  As  the  optical  power  grows  and  reaches  saturation,  the  figure 
shows  that  gain  and  mode  distortion  are  reduced.  Notice  that  at  about  n  =  200  passes, 
the  optical  mode  narrows  and  the  power  grows  rapidly  again,  finally  developing  a  more 
symmetric  Gaussian  mode  with  steady  power.  The  plot  of  |a(x,T)|  in  the  center  shows 
the  steady-state  optical  mode. 
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Figure  30.  Three-dimensional  simulation  results  for  the  TJNAF  10  kW  FEL, 
showing  optieal  mode  distortion  in  weak  fields,  and  the 
development  of  a  Gaussian  mode  at  saturation.  (From  [48].) 


2,  Current  NFS  Research:  Short  Rayleigh  Length  FELs 

In  the  previous  two  examples,  a  Rayleigh  length  of  Zg  =  0.4  was  ehosen  to 

maximize  the  overlap  between  the  eleetron  beam  and  the  optieal  mode.  In  a  high  average 
power  FEL,  it  may  be  desirable  to  operate  at  a  shorter  Rayleigh  length.  This  would  result 
in  a  smaller  mode  waist,  but  a  larger  spot  size  at  the  mirrors,  redueing  the  power  density 
on  the  mirrors.  Traditional  EEL  theory  implies  that  a  shorter  Rayleigh  length  would 
reduce  the  average  filling  factor,  and  hence  the  gain  and  steady-state  power,  but  mode 
distortion  can  counteract  this  effect.  We  have  used  our  three-dimensional  simulations  to 
study  short  Rayleigh  length  operation  for  the  proposed  10  kW  and  100  kW  upgrades  for 
the  TJNAE  EEL  at  Jefferson  Labs. 
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a.  TJNAFlOkWFEL 

Figure  31  shows  the  result  of  a  series  of  3dn  simulations  of  the  TJNAF 
10  kW  FEL.  The  weak-field  gain  and  steady-state  power  are  plotted  versus  the 
normalized  Rayleigh  length  Zq  .  The  design  value  of  Zg  =  0.4  is  indieated  by  the 
triangular  tiek  mark  at  the  bottom  of  the  graph.  The  weak-field  gain  remains  almost 
eonstant  as  the  Rayleigh  length  is  redueed  from  Zg  =  0.6  to  0.1.  The  gain  is  well  above 

threshold  {l/Q  »  15%)  for  the  entire  range  of  Zg  values  that  we  studied,  indieating  that  it 

may  be  possible  to  operate  at  an  even  shorter  Rayleigh  length.  As  the  Rayleigh  length  is 
redueed  from  the  design  value  of  Zg  =0.4  to  0.1,  the  steady-state  power  inereases  by 
50%.  Also  as  the  Rayleigh  length  is  redueed,  the  optimum  eleetron  phase  veloeity  Vg, 
indieated  at  each  point,  increases. 


0.1  0.2  0.3  0.4  0.5  0.6 

Normalized  Rayleigh  Length,  Zq 


Figure  31.  Three-dimensional  simulation  results  for  the  TJNAF  10  kW  FEL, 
showing  weak-field  gain  and  steady-state  power  versus  Rayleigh 
length.  As  the  Rayleigh  length  is  reduced,  the  weak-field  gain 
stays  fairly  constant,  and  the  steady-state  power  increases.  (Erom 
[48].) 
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b.  TJNAF 100  kWFEL 

The  design  for  the  proposed  100  kW  PEL  assumes  a  resonator  quality 
factor  Q  =  4.2  corresponding  to  approximately  21%  mirror  transmission,  and  edge  losses 
around  the  mirrors  of  1%  per  pass.  Even  with  such  large  losses,  the  intra-cavity  power 
may  be  too  much  for  the  mirrors  to  tolerate.  East  year  we  reported  the  results  of 
simulations  using  a  short  Rayleigh  length  resonator  to  reduce  the  power  the  density  on 
the  mirrors  [53]. 

Eigure  32  shows  the  shape  of  the  optical  mode  and  the  power  density  on 
the  mirrors  for  Rayleigh  length  Zg  =0.1  to  0.5.  The  horizontal  axis  corresponds  to  the 

longitudinal  axis  of  the  resonator,  normalized  to  the  undulator  length  L.  The  resonator 
length  is  1  \L  with  the  mirrors  located  at  t  =  -5  and  +6,  and  the  undulator  is  located 
between  r  =  0  and  1.  The  vertical  axis  is  the  optical  mode  radius  in  millimeters. 


Eigure  32.  Optical  mode  radius  versus  dimensionless  longitudinal  position  t 
for  an  optical  cavity  with  mirrors  located  at  r  =  -5  and  +6,  for 
various  Rayleigh  lengths.  As  the  Rayleigh  length  is  reduced,  the 
mode  radius  at  the  mirror  increases,  and  the  intensity  decreases. 
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Figure  32  shows  that  a  Rayleigh  length  of  Zg  =  0.3  gives  a  spot  size  at  the 
mirrors  of  1  cm,  corresponding  to  a  power  density  of  160  kW/cm^.  Reducing  the 
Rayleigh  length  from  Zg  =  0.3  to  0.1  increases  the  mirror  spot  size  to  approximately  1.7 
cm  and  decreases  the  power  density  on  the  mirrors  by  a  factor  of  3. 

A  short  Rayleigh  length  FEL  requires  a  short  undulator,  to  avoid  beam 
scraping  at  the  ends  of  the  undulator.  Also,  since  the  mode  waist  is  small  in  this  case,  a 
narrow  electron  beam  is  desirable.  Figure  33  shows  the  results  of  a  three-dimensional 
simulation  of  the  proposed  100  kW  FEL,  with  a  short  Rayleigh  length  given  by  Zg  =  0.1 , 

and  a  narrow  electron  beam  radius  given  by  =  0.1  in  the  x  and  y  directions.  In  this 

case  we  have  used  a  newer  simulation  called  3dbetan  that  includes  the  betatron  motion  of 
the  electrons.  The  betatron  oscillation  frequency  is  cOp  =1  over  the  undulator  length, 

with  the  electron  beam  focused  in  the  middle  of  the  undulator  at  Xp  =  0.5  .  The  beam’s 
angular  spread  Gg^^  =0.01  is  determined  by  the  matching  requirement  Gg^^  =C0pG^^. 
The  large  initial  electron  phase  velocity  Vg  =16  optimizes  the  extraction  efficiency  for 
the  short  Rayleigh  length  EEL.  The  dimensionless  mirror  radius  =  7.2  is  chosen  to 
allow  1%  edge  losses  per  pass  (e  =  0.0l),  and  the  dimensionless  radius  of  curvature 
=1.5  determines  the  Rayleigh  length. 

In  the  upper-left  of  Eigure  33,  a  slice  through  the  optical  wavefront 
|a(x,n)|  is  shown  as  it  evolves  over  n  =  200  passes.  After  about  50  passes,  the  optical 

field  begins  to  oscillate  between  two  different  modes  with  a  period  of  about  10  passes. 
The  optical  power  P{n)  and  gain  G[n)  in  the  lower  right  appear  to  oscillate  at  about  the 

same  frequency,  as  well  as  the  electron  phase  velocity  distribution  f(y,n)  in  the  lower- 
left.  The  final  wavefront  |a(x,y)|  in  the  upper-center  is  donut  shaped;  if  we  ran  the 
simulation  for  a  few  more  passes,  we  would  observe  a  different  mode  shape.  The  center 
plot,  |a(x,T  )| ,  presents  a  section  through  the  optical  wavefront  during  the  final  pass.  The 
superimposed  red  dots  show  the  sample  electrons  as  they  pass  through  the  undulator  from 
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T  =  0  to  1.  The  white  contour  line  shows  where  the  field  amplitude  drops  to  5%  of  its 
maximum  value,  as  indicated  on  the  color  scale  to  the  right.  The  complicated  pattern  of 
the  contour  line  confirms  the  presence  of  multiple  transverse  modes.  The  optical  wave 
amplitude  profile,  shown  in  yellow  at  the  mirrors  on  each  end,  indicates  that  power  is 
distributed  across  the  mirrors. 
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Figure  33.  Three-dimensional  simulation  results  for  the  proposed  100  kW 
TJNAF  FEL  with  a  short  Rayleigh  length,  Zg  =  0.1 ,  and  a  narrow 

electron  beam,  a  ^,=0.1,  showing  optical  mode  distortion  and 
limit-cycle  behavior. 


With  a  larger  electron  beam  radius,  multimodes  are  suppressed.  Figure  34 
shows  the  results  of  another  3dbetan  simulation  of  the  same  FEL,  with  a  longer  Rayleigh 
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length  given  by  Zq  =0.5  and  a  larger  electron  beam  radius  given  by  =0.5.  The 
dimensionless  current  density  has  been  reduced  to  7=14.4  in  order  to  keep  the  total 
current  jo p ^  constant.  The  electron  beam’s  angular  spread  Gg^^,  =0.25  has  been 
increased  to  satisfy  the  matching  condition  Gq^  „  =C0pG^,, .  With  the  longer  Rayleigh 
length,  the  optimum  initial  electron  phase  velocity  is  now  Vg  =  6.  The  dimensionless 
mirror  radius  that  corresponds  to  1%  edge  losses  is  now  =3.4,  and  the  dimensionless 
radius  of  curvature  is  =  1 .7  . 
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Figure  34.  Three-dimensional  simulation  results  for  the  proposed  100  kW 
TJNAF  FEL  with  a  long  Rayleigh  length,  Zg  =  0.5  ,  and  a  broad 

electron  beam,  g^,,  =0.5,  showing  the  development  of  a  stable 
Gaussian  mode. 
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In  Figure  34,  the  optieal  field  evolution  plots  |a(x,n)|  and  |a(x,T)|  show 

that  a  stable  Gaussian  mode  develops  over  n  =  32  passes.  Notiee  that  the  mirror  spot 
size  is  mueh  smaller  than  in  Figure  33,  and  the  power  is  eoncentrated  near  the  eenter  of 
the  mirror,  so  mirror  damage  is  more  likely  in  this  ease.  The  betatron  motion  of  the 
eleetrons  is  now  clearly  seen  in  the  red  dots  on  the  center-left  plot.  Some  of  the  sample 
electrons  make  excursions  far  outside  the  center  of  the  optical  mode,  where  they  won’t  be 
effectively  bunched.  As  a  result,  there  is  less  gain,  and  the  extraction  efficiency  is 
reduced  to  r]  «  0.8% . 

Figure  35  shows  a  plot  of  the  single-pass  extraction  efficiency  r]  versus 
the  electron  beam  radius  for  the  TJNAF  100  kW  FEL  parameters  with  a  short 
Rayleigh  length  Zq=0.1.  It  indicates  that  a  narrow  electron  beam  is  beneficial  for 
efficiency,  since  the  electrons  are  contained  well  inside  the  optical  mode. 
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Figure  35.  Three-dimensional  simulation  results  for  the  proposed  TJNAF  100 
kW  FEL  with  a  short  Rayleigh  length,  Zq  =  0.1 .  The  efficiency  r] 

is  plotted  versus  the  electron  beam  radius  The  extraction 

efficiency  increases  as  the  electron  beam  radius  is  reduced. 
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VI.  FOUR-DIMENSIONAL  SIMULATIONS 


The  simulations  described  in  the  last  two  chapters  follow  the  evolution  of  PEL 
optical  modes  in  either  the  longitudinal  or  transverse  directions.  They  do  not  require  too 
much  computer  processing  power  or  memory,  and  they  enable  us  to  isolate  and  study 
important  physical  effects  such  as  pulse  slippage,  coherence  evolution,  the  trapped- 
particle  instability,  diffraction,  mode  distortion,  and  wavefront  guiding,  as  discussed  in 
the  previous  chapters.  In  many  FELs,  either  longitudinal  or  transverse  mode  effects 
dominate,  depending  on  parameters  such  as  the  current  density,  the  electron  beam  radius, 
and  the  pulse  length.  But  in  general,  both  longitudinal  and  transverse  modes  are  present 
and  interact  with  each  other.  Eongitudinal  pulse  slippage  effects  tend  to  decrease  the 
gain,  while  transverse  mode  distortion  effects  tends  to  increase  the  gain.  It  is  sometimes 
necessary  to  study  both  types  of  effects  simultaneously.  For  that  reason,  we  have 
developed  four-dimensional  simulations  that  model  the  electron  pulse  and  the  optical 
wavefront  in  the  x,  y,  and  z  directions,  as  they  evolve  in  time.  Unlike  the  three- 
dimensional  simulations  described  in  the  last  chapter,  we  don’t  have  to  select  the 
resonance  parameter  v  in  advance;  the  four-dimensional  simulation  determines  the 
correct  value  of  v  self-consistently  through  mode  competition. 

Four-dimensional  simulations  may  be  useful  for  studying  optical  modes  in  a 
megawatt  class  FEE  with  a  short  Rayleigh  length.  However,  they  require  large  amounts 
of  memory  and  computer  time.  We  are  considering  porting  them  to  a  supercomputer,  or 
perhaps  building  a  special-purpose  computer  to  solve  this  problem.  We  are  also  studying 
new  techniques  for  reducing  the  computational  and  memory  requirements. 

A,  SHORT  PULSES 

For  modeling  short-pulse  FEFs  with  narrow  electron  beams,  we  have  developed  a 
multi-pass,  four-dimensional  simulation  called  4dpnf.  When  the  electron  beam  remains 
well  inside  the  optical  mode,  the  fdling  factor  is  small,  and  all 

electrons  interact  with  the  same  small  region  of  the  optical  field,  so  the  microscopic 
bunching  is  uniform  across  the  electron  beam.  The  distributions  of  electron  phases  ^ 
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and  phase  velocities  v  are  therefore  uniform  across  the  beam.  The  values  of  ^  and  v 
for  each  electron  are  independent  of  the  transverse  positions  x  and  y,  but  still  depend  on 
the  longitudinal  position  z  and  the  time  r  .  Instead  of  requiring  large  four-dimensional 
arrays  for  the  electron  microscopic  variables  C,  and  v  ,  only  two-dimensional  arrays  are 
needed.  This  reduces  the  major  memory  requirements  of  the  simulation.  Sample 
electrons  are  initially  distributed  uniformly  in  phase  from  0  to  2n  ,  but  may  have  a 
random  spread  in  phase  velocities  determined  by  beam  emittance  and  energy  spread. 

In  4dpnf,  a  complex-valued  three-dimensional  array,  a(x,y,z),  is  used  to  model 

the  optical  beam.  The  initial  optical  field  shape  is  parabolic  in  the  longitudinal  direction 
with  a  length  comparable  to  the  electron  pulse  plus  the  slippage  distance.  In  the 
transverse  direction,  the  initial  optical  field  amplitude  has  a  Gaussian  shape,  while  the 
optical  phase  is  curved  in  the  x-y  plane  so  that  natural  diffraction  focuses  the  beam  at 
position  along  the  undulator,  with  dimensionless  Rayleigh  length  Zq  =wl  and  mode 
waist  Wq  . 


A  three-dimensional  array  is  used  in  4dpnf  to  represent  the  electron  beam  current 
density  j{x,y,z') ,  which  is  assumed  to  have  a  parabolic  profde  in  the  x  andy  directions. 

Once  the  initial  electron  and  optical  beams  are  established,  the  electrons  evolve  in  phase 
space  according  to  the  pendulum  equation 


d\  (z-t,t) 
dx^ 


a(0,0,z,T)  cos[^  (z-T,T)  +  (t)  (0,0,z,t)], 


(52) 


while  the  optical  field  evolves  according  to  the  parabolic  wave  equation 


■-v^+A 
4  dx 


t{x,y,z,x)  =  -j{x,y,z-x){e 


(53) 


where  =8^  +  8^  with  x  and  y  again  normalized  to  the  characteristic  optical  mode 
radius  yjLX/n  ,  and  •  •)  indicates  an  average  over  all  electrons  at  each  optical  field  site. 

The  electron  pendulum  equation  is  evaluated  using  a  fourth-order  Runge-Kutta  method, 
while  the  parabolic  wave  equation  is  integrated  using  fast  Fourier  transforms  and 
diffraction  operators  [16]. 
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1. 


CEBAF  IR  FEL 


In  the  early  1990’s,  we  used  4dpnf  to  study  a  proposed  15  pm  IR  FEL  at  the 
Continuous  Eleetron  Beam  Accelerator  Eacility  (CEBAE,  now  known  as  the  Thomas 
Jefferson  National  Accelerator  Eacility  or  TJNAE)  [54],  This  author  was  deeply  involved 
in  planning  this  research,  developing  and  running  the  simulations,  and  interpreting  their 
output,  and  was  the  lead  author  of  a  paper  we  published  on  the  results  of  this  work  [55]. 

The  electron  beam  for  the  CEBAE  IR  experiment  had  an  energy  of  45  MeV,  with 
a  peak  current  of  72  A,  beam  radius  of  0.34  mm  and  pulse  length  of  0.54  mm.  The 
undulator  had  N  =  25  periods  of  length  kg  =  6  cm  each,  with  magnetic  field  strength 
B  =  4A  kG,  corresponding  to  undulator  parameter  A  =  1.76.  These  parameters  give  a 
peak  dimensionless  current  density  of  j  =  138,  dimensionless  beam  radius  of  =  0.125 
and  pulse  length  of  g^  =  1.4  .  A  dimensionless  Rayleigh  length  of  Zg  =  0.3  was  chosen  to 
maximize  the  theoretical  weak- field  gain,  with  the  mode  waist  at  =  0.5  .  The  resonator 
was  designed  to  have  2%  mirror  transmission,  corresponding  to  quality  factor  Q  =  50 , 
and  1%  edge  losses.  The  resonator  length  was  13  times  the  undulator  length,  so  a  mirror 
wide  enough  to  reflect  99%  of  the  light  must  have  a  radius  of  about  150  times  the 
electron  beam  radius.  That  disparity  in  sizes  requires  a  very  large,  fine  grid  in  the 
transverse  x-y  direction.  To  allow  for  a  more  reasonable  grid  size,  the  mirrors  were 
artificially  moved  closer  together,  in  this  case  to  dimensionless  separation  =  2  equal  to 

twice  the  undulator  length.  This  significantly  reduced  the  laser  spot  size  at  the  mirror. 
The  mirror’s  radius  of  curvature  was  increased  to  maintain  the  same  Rayleigh  length  over 
the  undulator.  If  the  mode  shape  is  the  same  over  the  EEL  interaction  region,  reducing 
the  length  of  the  free-space  propagation  region  outside  the  undulator  should  not  alter  the 
simulation  results. 

In  order  to  accurately  model  the  transverse  profile  of  the  electron  beam,  it  should 
exist  on  at  least  a  few  x-y  grid  points.  Eor  this  simulation,  we  increased  the  electron 

beam  radius  to  g^  =0.2  in  order  to  accommodate  a  reasonable  grid  size.  The  peak 
current  density  was  decreased  to  j  =  54  to  keep  the  same  total  current,  which  is 
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proportional  to  .  Since  the  eurrent  density  j  is  inversely  proportional  to  the  area  of 
the  eleetron  beam,  the  product  yaf  is  actually  independent  of  beam  size.  In  fact,  the 
PEL  single -mode  gain  is  proportional  to  y’F  oc  jo],  and  is  independent  of  the  eleetron 
beam  size.  If  the  transverse  mode  remains  undistorted,  then  the  value  of  has  no  effect 
on  the  gain.  However,  when  the  eleetron  beam  is  small,  transverse  mode  distortion  tends 
to  inerease  the  aetual  gain  above  the  single-mode  gain.  The  ehoiee  of  a^=0.2  is 
eonsidered  small  enough  to  allow  some  mode  distortion  if  it  should  oeeur. 

A  single  mode,  one-dimensional  simulation  in  t,  ineluding  the  filling  factor  F, 
predicts  about  50%  single-pass  gain  for  the  CEBAF  IR  EEE.  A  longitudinal  multimode, 
two-dimensional  simulation  in  (z,t)  shows  that  the  short  pulse  effeets  alone  reduee  the 
gain  to  about  25%.  With  a  transverse  multimode,  three-dimensional  simulation  in 
(x,y,t) ,  the  single-pass  gain  inereases  to  more  than  100%  beeause  of  mode  distortion.  A 

four-dimensional  simulation  in  (x,y,z,t),  whieh  ineludes  both  mode  distortion  and 

short-pulse  effeets,  gives  a  net  single-pass  gain  of  only  18%,  roughly  the  same  result 
given  by  the  two-dimensional  simulation.  The  results  of  the  4dpnf  simulation  over 
n  =  60  passes  are  shown  in  Figure  36.  The  output  is  similar  to  those  from  our  two- 
dimensional  simulations  described  in  the  ehapter  on  longitudinal  multimodes.  The  new 
intensity/eontour  plot  in  the  middle-right  is  the  evolution  of  a  cross-seetion  through  the 

optieal  wavefront  |a(x,n)|  at  y  =  z  =  0,  shown  at  the  end  of  eaeh  pass  through  the 

undulator.  The  final  wavefront  eross-seetion  after  n  =  60  passes,  shown  in  the  upper- 
right,  is  nearly  Gaussian.  In  the  middle  of  the  wavefront  is  a  dimple  that  is  about  the  size 
of  the  eleetron  beam,  shown  in  the  upper-right  as  superimposed  yellow  reetangle.  The 
dimple  is  eaused  by  the  interferenee  of  light  emitted  at  t  «  0.6  with  light  emitted  later  at 
t«0.9.  AtT«0.6,  where  eleetron  bunehing  begins,  the  new  light  generated  by  the 
narrow  eleetron  beam  of  radius  o^»0.2  has  a  Rayleigh  length  Zg^af  «0.04,  and 
rapidly  diffraets  outward.  The  phase  shift  of  the  generated  light  over  the  last  seetion  of 
the  undulator.  At  «  0.3 ,  ean  be  more  than  A(j)  «  -n  eausing  destructive  interferenee.  As 
a  result,  a  higher-order  mode  strueture  develops  over  many  passes.  The  final  mode  is  a 
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combination  of  the  fundamental  Gaussian  mode  and  a  higher-order  mode,  eausing  a 
dimple  on  the  scale  of  the  electron  beam.  This  dimple  disappears  when  the  electron  beam 
size  is  increased  beyond  ~  0.3 ,  beeause  the  phase  shift  is  insufficient  for  destructive 
interference  at  the  end  of  the  undulator. 
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Figure  36.  Four-dimensional  simulation  results  for  the  CEBAF  IR  FEL, 
showing  steady-state  gain  G{n)  in  weak  optical  fields,  and  the 

development  of  a  dimple  in  the  optical  wavefront,  |a(x,n)| .  (from 
[55].) 


The  theoretical  single-mode  gain  spectrum  for  these  parameters,  G  (v ) ,  is  shown 
in  the  lower  middle  of  figure  36  for  reference.  Above  the  gain  spectrum  is  the  optical 
power  spectrum  evolution,  P{y,n),  showing  the  initial  power  centered  at  resonance. 
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V  =  0.  After  n«15  passes,  the  laser  frequency  shifts  to  the  value  for  maximum  gain, 
corresponding  to  phase  velocity  v  «  2.6  .  The  final  power  spectrum  in  the  upper-center  is 
narrow,  but  is  Fourier- transform  limited  by  the  short  optical  pulse.  In  the  gain  evolution 
G{n)  shown  in  the  lower  right,  loss  initially  dominates  gain,  but  as  the  laser  frequency 

shifts  to  the  value  for  maximum  gain,  the  single-pass  gain  grows  to  about  18%. 
Considering  the  mirror  loss  and  the  effect  of  the  dimple,  this  is  in  rough  agreement  with 
the  weak-field  gain  of  25%  found  in  short-pulse  simulations  without  transverse  effects. 

Figure  37  shows  the  results  of  a  4dpnf  simulation  of  the  same  FEL,  allowed  to  run 
for  n  =  400  passes  to  reach  steady-state  saturation,  where  the  gain  each  pass  now  equals 
the  loss.  The  desynchronism  value  d  =  0.005  has  been  optimized  for  maximum  power  in 
strong  fields,  but  all  other  parameters  are  the  same  as  in  Figure  36.  The  optical  power 
spectrum  evolution  P{y  ,n)  in  the  center  of  Figure  37  shows  a  shift  to  phase  velocity 

V  «  271  for  optimum  gain  in  strong  fields,  as  we  saw  earlier  in  the  FEL  gain  surfaces. 
The  final  optical  mode  shape  |a(x,n)|  has  no  dimple,  because  the  strong  fields  reduce  the 
phase  shift  that  is  needed  for  interference. 
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Figure  37.  Four-dimensional  simulation  results  for  the  CEBAF  IR  FEL, 
showing  steady-state  optical  power  P{n)  and  the  development  of 

a  smooth  Gaussian  wavefront  |a(x,n)|.  (From  [55].) 


B.  LONG  PULSES 

For  modeling  long-pulse  FELs,  we  have  developed  a  multi-pass,  four-dimensional 
simulation  called  4dwrn,  which  includes  the  betatron  motion  of  the  electrons  and  does  not 
assume  a  narrow  electron  beam.  A  complex -valued  three-dimensional  array,  a(x,y,z), 
is  used  to  model  the  optical  beam  inside  a  window  of  width  a  few  slippage  distances  or 
less.  A  three-dimensional  array  represents  the  electron  beam  current  density  7(x,y,z), 

which  is  assumed  to  have  a  parabolic  profile  in  the  x  andy  directions.  One-dimensional 
arrays  are  used  to  store  the  microscopic  phase  and  phase  velocity  of  each  sample 
electron,  and  two-dimensional  arrays  represent  their  x  and  y  positions  as  they  execute 
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betatron  oscillations.  As  before,  the  pendulum  equation  (52)  is  used  to  follow  the 
microscopic  motion  of  the  electrons  in  phase  space,  and  the  parabolic  wave  equation  (53) 
describes  the  evolution  of  the  light. 

1.  Boeing  1  kW  FEL  Oscillator 

In  1996,  we  used  4dwrn  to  study  the  Boeing  1  kW  FEL  oscillator,  discussed 
earlier  in  the  chapter  on  longitudinal  modes.  This  author  played  a  key  role  in  planning 
this  research,  developing  and  running  the  simulations,  and  interpreting  their  output,  and 
was  the  lead  author  of  a  paper  we  published  on  the  results  of  this  work  [56].  Figure  38 
shows  the  output  from  a  4dwrn  simulation  of  the  Boeing  FEL  oscillator  over  n  =  50 
passes.  The  dimensionless  current  density  is  7=1800,  with  an  electron  beam  radius 
Gg  =  0.5  .  The  electrons  have  a  Gaussian  energy  spread  given  by  g^  =  1 1  and  an  angular 
spread  Gg  =1.4.  The  electrons  have  a  betatron  frequency  cOp  =  8.4 ,  or  slightly  more  than 
one  betatron  oscillation  along  the  N  =  220  period  undulator.  The  dimensionless 
Rayleigh  length  is  Zg  =  0.5  ,  with  «20%  mirror  outcoupling  [Q  =  5)  and  1%  edge  losses 

(e  =  0.0l). 

The  power  evolution  P{n)  in  the  lower  right  of  Figure  38  shows  that  large 

optical  fields  develop  rapidly  inside  the  high-gain  FEL  (7»l),  creating  deep  potential 

wells  in  phase  space.  The  bunched  electrons  undergo  synchrotron  oscillations,  and  the 
trapped-particle  instability  leads  to  sideband  formation.  The  sidebands  are  amplified 

each  pass,  modulating  the  field  envelope  |a(z,n)|  shown  on  the  left,  and  producing  a 

broad,  chaotic  optical  power  spectrum,  shown  in  the  center  plot,  P{y,n).  The  final 

spectral  width  is  Ak/A,  «  Av  jlnN  ~  3.5%  .  The  simulation  predicts  a  steady-state  output 

power  of  F  «  8x  10"^  in  our  dimensionless  units,  which  corresponds  to  about  250  W,  after 
accounting  for  the  overall  system  duty  factor. 
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Figure  38.  Four-dimensional  simulation  results  for  the  Boeing  1  kW  FEL 
oseillator,  showing  the  development  of  a  broad,  ehaotie  optieal 
power  spectrum  P(y  ,n).  (From  [56].) 


One  way  to  prevent  the  trapped-particle  instability  is  to  taper  the  undulator.  This 
should  cause  about  half  of  the  electrons  to  become  untrapped,  which  would  reduce  the 
synchrotron  oscillations.  Figure  39  shows  another  4dwrn  simulation  of  the  Boeing  FEL 
oscillator,  with  a  10%  undulator  taper,  corresponding  to  phase  acceleration  5  =8871:  . 
The  theoretical  weak-field  gain  spectrum  G  (v )  in  the  lower  center  has  numerous  spikes 
because  of  the  large  undulator  taper. 
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Figure  39.  Four-dimensional  simulation  results  for  the  Boeing  1  kW  FEL 
oseillator  with  a  tapered  undulator,  5  =8871:  ,  showing  the 
suppression  of  sidebands  and  the  development  of  a  narrow  optical 
power  spectrum  P{y  ,n).  (From  [56].) 


In  Figure  39,  we  see  that  the  sideband  growth  is  reduced  significantly,  and  the 
spectral  width  is  now  only  Ak/A,  »  1% .  The  average  output  power  is  increased  to  430  W, 

but  the  power  evolution  P{n)  shows  significant  oscillations.  This  may  be  due  to 

sporadic  attempts  at  sideband  formation  and  interference  of  sideband  frequencies. 
Increasing  the  outcoupling  or  the  system  duty  factor  would  be  necessary  to  achieve  1  kW 
average  power. 
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C.  CURRENT  NFS  RESEARCH:  MEGAWATT  EEL 

At  the  Naval  Postgraduate  School,  we  are  working  on  the  design  of  a  compact, 
high  average  power  PEL  for  ship  defense  in  collaboration  with  Jefferson  Labs  and 
Advanced  Energy  Systems  [57].  The  laser  must  produce  at  least  1  MW  of  average  power 
at  a  wavelength  near  1  pm  for  propagation  through  the  marine  atmosphere  and 
destruction  of  incoming  missiles.  With  25%  outcoupling,  the  intra-cavity  power  would 
be  4  MW.  The  entire  system  is  intended  to  fit  inside  a  box  this  is  approximately 
12m  X  4m  X  4m ,  so  the  resonator  can  only  be  12  m  long.  Eor  a  typical  Rayleigh  length  of 
Zg  =  5  m,  the  laser  spot  size  on  the  mirrors  would  be  only  a  few  mm  diameter.  This  spot 

size  would  give  a  power  density  on  the  mirrors  of  about  30  MW/cm  ,  which  would 
destroy  the  mirrors.  We  have  proposed  an  optical  resonator  with  a  very  short  Rayleigh 
length  of  Zq=1.8  cm  to  increase  the  mirror  spot  size  and  reduce  the  intensity  to 
210  kW/cm  .  Cooled  sapphire  mirrors  may  be  able  to  withstand  this  intensity  without 
damage. 

The  design  calls  for  a  short  A  =  20  period,  T  =  60  cm  length  undulator,  to 
enhance  the  extraction  efficiency  and  avoid  beam  scraping  at  the  ends  of  the  undulator. 
Eor  these  parameters,  the  dimensionless  Rayleigh  length  is  Zg  =Zo/E  =0.03,  compared 
to  a  typical  EEL  with  a  Rayleigh  length  of  Zg  =  0.3 .  The  rms  undulator  parameter  for  this 
design  is  K  =  2. 

The  optical  mode  radius  at  its  waist  in  the  middle  of  the  proposed  resonator  is 
only  Wg  =  0.08  mm.  With  such  a  short  Rayleigh  length,  the  beam  will  expand  rapidly,  to 

a  spot  size  of  2.5  cm  at  the  mirrors,  approximately  300  times  greater  than  at  the  waist. 
This  size  disparity  implies  that  a  grid  of  at  least  4096x4096  elements  in  the  transverse 
(x,y)  plane  is  required  to  represent  the  optical  beam  throughout  the  resonator.  The  grid 

dimensions  are  determined  by  the  requirement  that  it  must  cover  the  entire  optical 
wavefront  at  the  mirrors,  yet  still  represent  the  much  smaller  wavefront  at  the  waist  with 
>  16  grid  elements.  Note  that  the  EETs  used  in  our  propagation  algorithm  require  that 
the  grid  dimensions  be  integral  powers  of  2. 
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The  superconducting  accelerator  will  provide  an  electron  beam  with  an  energy  of 
140  MeV  and  an  average  current  of  0.6  A.  The  dimensionless  current  density  is  j  =  200 , 
so  this  is  a  high-gain  PEL.  The  electron  beam  will  be  recirculated  to  enhance  overall 
system  efficiency  and  reduce  shielding  requirements  for  the  beam  dump.  An  electron 
beam  radius  of  =  0.08  mm  would  optimize  the  overlap  with  the  optical  mode  in  the 
center  of  the  undulator.  This  corresponds  to  a  normalized  electron  beam  radius  of 
^  =  0.2 ,  so  there  could  be  significant  optical  mode  distortion. 

The  electron  pulse  length  is  4  =0.3  mm,  or  a  normalized  length  of  g^=15 

slippage  distances,  so  pulse  slippage  effects  should  be  minor.  However,  longitudinal 
modes  should  be  included  to  study  coherence  evolution,  and  to  allow  the  simulation  to 
self-consistently  determine  the  laser  operating  frequency.  To  accurately  model  the 
longitudinal  pulse  for  a  reasonable  number  of  time  steps  through  the  undulator,  at  least 
1024  sites  are  required  in  the  z  direction.  With  the  requirement  for  4096  x  4096  elements 
in  the  (x,y)  plane  at  each  of  the  longitudinal  sites,  assuming  8  bytes  for  each  double¬ 
precision  floating-point  number,  a  256  GB  array  is  needed  to  represent  the  amplitude  and 
phase  of  the  optical  field  at  every  point  inside  the  optical  cavity.  This  is  much  greater 
than  the  capacity  of  currently  available  personal  computers  and  desktop  workstations. 

To  avoid  this  memory  limitation,  we  could  artificially  move  the  mirrors  in  closer, 
as  we  have  done  in  some  of  three-dimensional  simulations  described  earlier.  Another 
possibility  is  to  use  a  re-gridding  technique.  At  the  end  of  the  undulator,  the  mode  radius 
is  only  =  1.3  mm,  about  16  times  the  waist  radius.  So  inside  the  undulator  we  could 
use  a  grid  of  only  256x256  elements  in  the  (x,y)  plane,  requiring  only  1  GB.  Outside 

the  undulator,  when  the  optical  beam  expands  rapidly  and  no  longer  interacts  with  the 
narrow  electron  beam,  we  could  map  the  optical  wavefront  onto  a  coarser  grid,  the  same 
size  (256x256)  but  with  each  element  representing  a  larger  area.  We  would  use  this 

coarse  grid  to  propagate  the  wavefront  out  to  the  mirror  and  back,  and  then  map  it  back 
onto  the  fine  grid  when  it  re-enters  the  undulator.  One  issue  we  are  studying  is  how  to 
map  the  wavefront  from  one  grid  to  another  without  introducing  numerical  artifacts.  We 
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believe  re-gridding  won’t  affect  the  fundamental  mode,  but  might  cause  distortion  of 
higher-order  modes. 
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VII.  CONCLUSION 


The  results  presented  in  this  dissertation  cover  10  years  of  PEL  research  at  the 
Naval  Postgraduate  School.  This  work  has  added  to  the  understanding  of  important 
effects  in  FELs,  such  as  the  trapped-particle  instability,  pulse  lethargy,  coherence 
evolution,  limit-cycle  behavior,  gain  degradation  due  to  electron  beam  quality,  optical 
mode  distortion,  and  optical  guiding.  It  has  helped  to  explain  new  physical  effects  such 
as  electron  beam  energy  modulation.  It  has  guided  the  development  of  new  FEE 
concepts  such  as  multi-stage  optical  klystrons,  inverse-tapered  and  step-tapered 
undulators,  and  short  Rayleigh  length  resonators.  It  has  played  an  important  role  in  the 
design  of  many  existing  and  proposed  FEE  experiments,  including  the  Boeing  1  kW 
oscillator  and  APEE  amplifier,  the  Stanford  FIREFEY  and  SCA  FEEs,  the  SEEENE  FEE 
klystron,  the  SLAC  X-ray  FEE,  the  CEBAF  IR  FEE,  and  the  TJNAF  FEEs. 

Our  most  recent  work  has  concentrated  on  the  development  of  high  average- 
power  FEEs  for  naval  applications.  We  have  worked  closely  with  colleagues  at  TJNAF 
as  they  upgrade  their  present  FEE  from  2  kW  to  10  kW,  and  eventually  to  100  kW 
average  power.  We  have  visited  their  laboratory  numerous  times,  attended  workshops 
and  upgrade  review  panels,  and  produced  several  joint  publications  and  conference 
presentations.  An  important  research  goal  has  been  to  enhance  their  FEE  extraction 
efficiency  while  minimizing  the  induced  energy  spread  of  the  electrons,  to  enable  beam 
recirculation.  To  reach  this  goal,  they  are  considering  new  undulator  designs  such  as 
inverse-tapered  and  step-tapered  undulators.  The  research  results  presented  in  Chapter 
IV  of  this  dissertation  have  been  useful  in  evaluating  these  new  undulator  designs.  Our 
simulations  showed  that  a  small  negative  linear  or  step-taper  produces  the  best  efficiency, 
but  only  a  slight  improvement  over  no  taper,  in  contrast  to  previously  published  results. 

A  high  average-power  FEE  may  require  a  short  Rayleigh  length  resonator  to 
avoid  mirror  damage.  Classical  FEE  theory  says  that  the  Rayleigh  length  should  be  about 
one  third  of  the  undulator  length  to  optimize  the  overlap  between  the  electron  beam  and 
the  optical  mode.  However,  the  research  results  presented  in  Chapter  V  of  this 
dissertation  show  that  it  may  be  possible  to  operate  an  FEE  at  a  much  shorter  Rayleigh 
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length,  one  tenth  of  the  undulator  length  or  less,  without  a  significant  reduction  in  gain. 
In  fact,  our  simulation  results  shown  in  Figure  31  predict  that  the  steady-state  power 
actually  increases  as  the  Rayleigh  length  is  reduced.  We  believe  this  effect  is  due  to 
optical  mode  distortion,  as  seen  for  example  in  Figure  33. 

At  the  Naval  Postgraduate  School,  we  are  currently  designing  a  compact, 
megawatt  average-power  FEL  for  ship  defense,  in  collaboration  with  TJNAF  and 
Advanced  Energy  Systems.  To  satisfy  the  size  requirements  without  damaging  the 
mirrors,  the  optical  cavity  may  need  to  have  a  very  short  Rayleigh  length.  In  that  case, 
our  present  3D  and  4D  simulations  will  need  to  use  very  large  transverse  grids, 
overwhelming  the  memory  capacity  of  desktop  computers.  To  avoid  this  memory 
limitation,  we  are  developing  a  re-gridding  technique,  discussed  in  Chapter  VI. 

The  original  research  by  the  author  presented  in  this  dissertation  has  played  an 
important  role  in  the  development  of  EEEs  over  the  past  10  years.  Using  numerical 
simulations,  this  work  has  helped  to  bridge  the  gap  between  EEE  theory  and  experiments, 
to  understand  important  physical  effects  in  EEEs,  and  to  aid  in  the  design  of  new  EEEs. 
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